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(54) Method of detection of angular velocity and torque In an Internal comlwstlon engine 



(57) An engine in which an elapsed time Ta(i) of 30** 
crank angle near top dead center of the compression 
stroke and an elapsed time Tb(i) of 30* crank angle near 
90** after top dead center of the compression stroke are 
found. When the elapsed time Ta(i) increases due to the 
torsional vibration of the engine drive system, the 
detected elapsed time Tb(i) is conrected to reduce it. 
while when the elapsed time Ta(i) decreases, the 
detected elapsed time Tb(i) is corrected to increase it 
The drive force or torque generated at each cylinder is 
calculated fron the thus downward corrected or upward 
corrected elapsed time Tb(i) and the elapsed time Ta(i). 
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Description 

BACKGROUND OF THE INVENTION 
5 1 . Reld of the Invention 

The present invention relates to a method of detection in an internal combustion engine. 
2. Desaiption of the Related Art 

10 

Known in the art is an internal combustion engine which finds a first angular velocity of the crantehaft in the time 
required for the crankshaft to rotate from 30*" to 60*" after top dead center of the compression stroke from this period, 
finds a second angular velodt/ of the crankshaft in the time required for the cranksliaft to rotate from 90*" to 120"* after 
top dead center of the compression stroke from this time, finds the torque generated by a cylinder from the square of 
IS the first angular velocity and the square of the second angular velocity, and calculates the amount of fluctuation of the 
torque from the amount of fluctuation of the generated torque (see Japanese Examined Patent Publication (Kokoku) No. 
7-33809). 

That is. when combustion is performed in a cylinder, the combustion pressure causes the angulcu: velocity of the 
crankshaft to rise from a first angular velocity coa to a second angular velodty cob. At this time, if the moment of inertia 

20 of rotation of the engine is I. the corrbustion pressure causes the kinetic energv to rise from (1/2) - Icoa^ to (1/2) • loob^. 
Roughly speaking, the amount of rise of tiie kinetic energy (1/2) • I • (a>b ^- cf>a ) causes a torque to be generated, so 
the generated torque becomes proportional to (o)b^-coa^). Therefore, the generated torque is found from the difference 
between the square of the first angular velocity oca and the square of the second angular velocity cob and. therefore, in 
the above-mentioned internal combustion engine, the amount of fluctuation of tiie torque is calculated from the thus 

25 found generated torque. 

However, if the generated torque is calculated based on the angular velocities oca and cob in this way. when for 
example the engine drive system experiences torsional vibration, the generated torque calculated based on the angular 
velocities oa and cob will no longer express the true generated torque. That is. when the engine drive system does not 
experience a torsional vibration, the second angular velocity cob increases from the first angular velocity coa by exactiy 

30 the amount of increase of the angular velocity caused by the combustion pressure. As opposed to this, when the eigine 
drive system experiences a torsional vibration, the second angular velcx:ity cob will include in addition to the amount of 
increase of tiie angular velocity caused by the combustion pressure the amount of change of the angular velocity 
caused by the torsional vilsration of the engine drive system in the period from detection of the first angular velocity coa, 
to detection of the second angular velocity cob. For example, rf the angular velocity increased due to the torsional vibra- 

35 tion of tiie engine drive system in the period from detection of the first angular velcx:ity coa to detection of tiie second 
angular velocity db, the amount of increase of the second angular velocity cob with respect to the first angular velocity 
coa will include in addition to the amount of increase of the angular velocity due to the combustion pressure the amount 
of increase of the angular velocity due to the torsional vibration of the engine drive system. Therefore, in this case, in 
so far as the amount of increase of the angular velocity due to the torsional vibration of the engine drive system is not 

40 subtracted from the second angular velocity cob, the difference between the square of the first angular velocity coa and 
the square of the second angular velocity cob will not express the generated torque. 

However, in the above-mentioned internal combustion engine, no consideration at all is given to the amount of 
change of the angular velocity due to the torsional vibration of tiie engine drive system and therefore, in the alx>ve-men- 
tioned internal combustion engine, when the engine drive system experienced torsional vibration, there was the prob- 

45 lem tiiat it was not possible to detect the true generated torque. 

SUMMARY OF THE INVENTION 

An object of the present invention is to provide a method of detection in an engine capable of accurately detecting 
50 an operating parameter of an engine such as a driving force of an engine, an output torque of an engine, and an amount 
of fluctuation of the output torque of an engine. 

According to tiie present invention, there is provKled a method of detection in an internal combustion engine, com- 
prising the steps of setting a first crank angle range in a crank angle region from the end of a compression stroke to the 
beginning of an expansion stroke, setting a second crank angle range in a crank angle region in the middle of the expan- 
55 ston stroke a predetermined crank angle away from the first crank angle range, detecting a first angular velocity of the 
crankshaft in the first crank angle range, detecting a second angular velocity of tiie aankshaft in the second crank angle 
range, finding the amount of change of the angular velocity between cylinders from the difference of the first angular 
velocity of a cylinder previously performing comt>ustion and the first angular velocity of a cylinder next performing com- 
txjstion, correcting the second angular velocity of the cylinder previously peribrming the combustion in the downward 
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direction when the amount of change of the angular velocity between cylinders has increased, correcting the second 
angular velocity of the cylinder previously pertornfiing the combustion in the upward direction when the amount of 
change of the angular velocity between cylinders has decreased, and finding the drive force generated from a corre- 
sponding cylinder based on the first angular velocity and the conrected second angular velocity. 

5 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention nnay be more fully understood from the description of the preferred embodiments of the 
invention set forth below together with the accompanying drawings, in which: 

10 

Fig. 1 is an overall view of an internal comtxjstion engine; 

Rg. 2 is a view of a map of the basic fuel irijection time; 

Rg. 3 is a view of the amount of generation of NOx €urid torque fluctuation; 

Rg. 4 is a view of a map of a lean correction coefficient; 
IS Rg. 5 is a view of a map of a lean limit feedback correction coefficient; 

Rgs. 6A and 6B are time charts of the changes in the elapsed times Ta(i) and Tb(i) of 30^ crank angle; 

Fig. 7 is a time chart of the changes in the elapsed time Ta(i) of 30'' crank angle; 

Fig. 8 is a time chart of the changes in the elapsed times Ta(i) and Tb(i) of 30'' crank angle; 

Fig. 9 is a time chart of the changes in the elapsed times Ta(0 and Tb(i) of 20"* crank angle; 
20 Rg. 10 is a time chart of the changes in the elapsed time Ta(0 of 30** crank angle; 

Rg. 11 is a flowchart of the interruption routine; 

Rg. 12 is a flowchart for calculating the elapsed times Ta(i) and Tb(i); 

Figs. 13 to 1 5 are flowcharts for checking the permission for calculation of the torque; 

Fig. 16 is a time chart of the changes of tiie elapsed time Ta(i) and the changes of the flags XMXREG and XMN- 
25 REG; 

Fig. 1 7 is a flowchart for cateulating the torque; 

Rgs. 1 8 and 1 9 are flowcharts for calculating the ratios KTa(i) and KTt)(i); 
Rg. 20 is a flowchart for processing of tiie counter GDLNIX; 
Rg. 21 IS a view of the timings for calculation of various values; 
30 Figs. 22 A and 228 are views of a target torque fluctuation value; 

Figs. 23A and 23B are views of tiie of fluctuation amount judgement values DH(n) arid DL(n) and the levels of ^ 
torque fluctuation LVLH(n) and LVLL(n) ; 
Rg. 24 is a flowchart showing a main routine; 

Rgs. 25 and 26 are flowcharts for calculating the torque fluctuation value; 
35 Rg. 27 Is a flowchart of the calculation of a lean limit feedback correction coefficient; 
Rg. 28 is a flowchart for calculating the fuel injection time; 

Rg. 29 is a view of the relationshf) of the anplitude of the torsional vibration of the crankshaft and the engine speed 
N; 

Rgs. 30A and 30B are views of changes in the angular velocity; 
40 Rg. 31 is a view of the relationship of the aankshaft position and the amplitude of the torsional speed; 

Rg. 32 is a flowchart of an interruption routine; 

Fig. 33 is a flowchart for calculating the lapsed times Ta(i) and Tb(0; 

Rgs. 34 and 35 are flowcharts for calculating the torque; 

Fig. 36 is a flowchart showing a main routine; 
45 Figs. 37 and 3d are flowcharts for calculating the torque fluctuation value; 

Rg. 39 is a flowchart for calculating the fuel injection time; 

Rg. 40 is a time chart of the changes in the elapsed times Ta(t) and Tb(i} of 30** crank angle; 
Figs. 41 A and 41 B are time charts of the changes in the elapsed times Ta(i) and Tb(i) of 30"* crank angle; 
Rg. 42 is a view of the amount of increase of the elapsed time; 
so Fig. 43 is an enlarged side view of a rotor; 

Rg. 44 is a view of the amount of increase of the elapsed time; 
Fig. 45 is a flowchart of an interruption routine; 

Rg. 46 is a flowchart for checking the permission for calculation of the torque; 
Rgs. 47 to 49 are flowcharts for calculating the torque; 
55 Rgs. 50 and 51 are flowcharts for calculating the ratios KTa(i) and KTb(i); 
Fig. 52 is a flowchart for processing of the counter GDLNIX; 
Fig. 53 is a view of the timings for calculation of various values; 
Fig. 54 is a flowchart showing a main routine; 

Rgs. 55 and 56 are flowcharts for calculating the torque fluctuation value; 
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Rg. 57 is a flowchart for calculating the fuel injection time; 

Fig, 58 is a time chart off the changes in the elapsed time Ta(i) of 30* crank angle and the elapsed time Tb(i) of 50" 
crank angle; 

Fig. 59 is a time chart of the changes in the elapsed time Ta(i) of 30"* crank angle and the elapsed time Tb(i) of 50** 
5 crank angle: 

Fig, 60 is a time chart of the changes in the elapsed time Ta(Q of SO"" crank angle and the elapsed time Tb(i) of 
crank angle; 

Fig. 61 is a flowchart of an interruption routine; 
Rg. 62 is a flowchart for calculating the elapsed times Ta(i) and Tb(i); 
10 Rgs. 63 to 65 are flowcharts fbr checking the permission for calculation of the torque; 
Rg. 66 is a flowchart for calculating the torque; and 
Rg. 67 is a flowchart showing a main routine, 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

IS 

Referring to Rg. 1 . 1 shows an engine body provided with four cylinders consisting of the No. 1 cylinder #1 , No. 2 
cylinder #2, Na 3 cylinder #3, and No. 4 cylinder #4. The cylinders #1» #2. #3. and #4 are respectively connected 
through the corresponding intake pipes 2 to a surge tank 3. In the intake tubes 2 are provided fuel injectprs 4 fbr inject- 
ing fuel toward the corresponding intake ports. The surge tank 3 is connected through an intake duct 5 to an air cleaner 

20 6. In the intake duct 5. a throttle valve 7 is arranged. On the other hand, the cylinders #1 . 1f2, #3. and #4 are connected 
through an intake manifold 8 and an exhaust pipe 9 to a casing 11 aooomnrKxiating an NOx absort)ent 10. This NOx 
absorbent 10 has the function of absorbing the NOx included in the exhaust gas when the air-fuel ratio is lean and dis^ 
charging the atisorbed NOx when the air-fuel ratio is the stoichiometric air-fuel ratio or rich. 

The electronic control unit 20 is comprised of a digital computer and is provided with a ROM (read only memory) 

25 22. a RAM (random access memory) 23, a CPU (microprocessor) 24, a backup RAM 25 connected to a constant power 
supply, an input port 26. and an output port 27 connected witii each other by a bidirectional bus 21 . The output shaft of 
the engine 12 has attached to it a rotor 13 with outer teeth. A crank angle sensor 14 comprising an electronr»gnetic 
pickup is arranged facing the outer teeth of the rotor 13, As shown in Rg. 1, in this ennbodiment. the rotor 13 has an 
outer tooth formed on its periphery at every 30*" crank angle and. fbr example, has part of the outer teetii removed for 

30 detecting the top dead center of the compression stroke of the No. 1 cylinder. Therefore, except for the portion where 
the outer teeth are removed, that is, the non-tooth portion, the crank angle sensor 1 4 generates an output pulse every 
. time the output shaft 12 turns by 30" crank angle. This output pulse is input to the input port 26. 

The surge tank 3 has attached to it a pressure sensor 1 5 for generating ah output voltage proportional to the abso-. 
lute pressure in the surge tank 3. The output voltage of this pressure sensor 15 is input through a conresponding AD 

35 converter 28 to the input port 26. Further, the throttie valve 7 has attached to It an kJle switch 16 fbr detecting when the 
throttle valve 7 is at the idling opening position. The output signal of this idle switch 16 is input to the input port 26. Fur- 
ther, the intake manifoki 8 has disposed in it an air-fuel ratio sensor (O2 sensor) 1 7 for detecting the air-fuel ratio. The 
output signal of this air-fuel ratio sensor 1 7 Is input through the corresponding AD converter 28 to tiie input port 26. On 
the otiier hand, the output port 27 is connected tiirough the corresponding drive circuit 29 to the fuel injectors 4. 

40 In the internal combustion engine shown in Rg. 1 , the fuel injection time TAU is calculated based on the following 
equation: 

TAU = TP • FLEAN • FLLFB • FAF+TAUV 

45 Here. TP shows a basic fuel injection time. FLEAN a lean correction coefficient. FLLFB a lean limit feedback cor- 
rection coefficient, FAF a stoichfometric air-fuel ratio feedback correction coefficient, and TAUV an invalid injection time. 

The basic fuel injection time TP shows the injection time required fbr making the air-fuel ratio the stoichiometric air- 
fuel ratio. This bask: fuel injection time TP is foiond from experiments. This ksasic fuel Injection time TP is stored in the 
ROM 22 in advance in the form of a map shown in Fig. 2 as a function of the absolute pressure PM in the surge tank 3 
so and the engine speed N. 

The lean correction coeff iderrt FLEAN is a con-ection coefficient for making the air-fuel ratio a target lean air-fuel 
ratio. This lean conrection coefficient FLEAN is stored in advance in tiie ROM 22 in the form of the map shown in Rg. 4 
as a function of the absolute pressure PM in the surge tank 3 and the engine speed N. 

The lean limit feedback correction coefficient FLLFB is a connection coefficient for maintaining the air-fuel ratio at 
55 the lean limit In this embodiment according to the present invention, the learning region for the lean air-fuel ratio feed- 
back control for the absolute pressure PM in the surge tank 3 and the engine speed N is divided into nine regions as 
shown in Rg. 5 for example. Lean limit feedback conrection coefficients FLLFB 1 1 to FLLFB 33 are set for the learning 
regions. 

The stoichiometric air-fuel ratio feedback correction coefficient FAF is a coefficient for maintaining the air-fuel ratio 
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at the stoichiometric air-fuei ratio. The stoichiometric air-fuel ratio feedback correction coefficient FAF is controlled 
based on the output signal of the air-fuel ratio sensor 1 7 so as to maintain the air-fuel ratio at tiie stoichiometric air-fuel 
ratio. At this time, the stoichiometric air-fuel ratio feedback correction coefficient FAF varies sut^stantially about 1 .0. 

The lean correction coefficient FLEAN is set in accordance with the operating state of the engine for the operating 
5 region enclosed by the broken lines in Rg. 4. in this operating region, the air-fuel ratio is maintained at the target lean 
air-fuel ratio. As opposed to this, in the operating region outside the region enclosed by the broken line in Rg. 4. the air- 
fuel ratio is maintained at tiie stoichiometric air-fuel ratia When the air-fuel ratio is to be maintained at the stoichiometric 
air-fuel ratk>, the lean correction coefficient FLEAN and the lean linm't feedback correction coeff ictent FLLFB are fixed at 
1 .0 and the stoichiometric air-fuel ratio feedback correction coefficient FAF is controlled based on the output signal of 
10 the air-fuel ratio sensor 1 7. 

On the other hand, when the air-fuel ratio is to be maintained at the target lean air-fuel ratio, the stoichiometric air- 
fuel ratio feedback correction coefficient FAF is fixed at 1 .0. that is. the feedback control based on the output signal of 
the air-fuel ratio sensor 17 is stopped, and tiie lean correction coefficient FLEAN and tiie lean limit feedback correction 
coefficient FLLFB are used to control tiie air-fuel ratio to tiie target lean air-fuel ratia 
75 Next, the lean limit feedback control will be explained refem'ng to Rg. 3. Rg. 3 shows the relationship between the 
Emiount of fluctuation of the torque of the engine output and the amount of generation of NOx and the air-fuel ratio. The 
leaner the air-fuel ratio, the smaller the fuel consumption rate. Further, the leaner the air-fuel ratio, tiie smaller the 
amount of generation of NOx. Therefore, viewed from these points, the air-fuel ratio should desirat)ly be made as lean 
as possible. Note, however, that when tiie air-fuel ratio becomes leaner than a certain extent, tiie combustion becomes 
20 unstable and, as a result, as shown in Fig. 3, the amount of fluctuation of the torque becomes large. Therefore, in this 
embodlmerrt according to the present invention, as shown in Rg. 3, the air-fuel ratio ismaintained in the air-fuel ratio 
control region where the torque fluctuation starts to increase. 

That is, explaining this in more detail, the lean correction coefficient FLEAN is determined so that the air-fuel ratio 
becomes the middle of tiie air-fuel ratio control region shown in Rg. 3 when the lean limit feedback correction coefficient 
25 FLLFB is made 1 .0. On the other hand, the lean limit feedttack correction coeffident FLLFB is controlled to within tiie 
torque fluctuation control region shown in Rg. 3 in accordance with the amount of fluctuation of the torque. When the 
amount of fluctuation of the torque becomes larger, the lean limit feedback correction coeff ident FLLFB is increased, 
that is. the air-fuel ratio is made smaller, while when the amount of fluctuation of the torque becomes smaller, the lean 
limit feedback correction coeffident FIJ.FB is reduced, tiiat is, the air-fuel ratio is made larger. In this way, the air-fuel 
30 ratio is controlled to within the air-fuel ratio control region shown in Rg. 3. / 

Note that, as will be understood from a comparison of Rg. 4 and Rg. 5. the lean limit feedback correction coefficient 
FLLPB is set to substantially the same region as the engine operating region where the leEih. correction coefficient 
FLEANisset 

\i the amount of fluctuation of tiie torque is controlled to wiUiin the torque fluctuation control region shown in Rg. 3. 

35 an excellent drivability of the vehide may be ensured while the fuel consumption rate arxj the amount of generation of 
NOx can be greatiy reduced. However, since the amount of fluctuation of the torque is controlled to be witiiin the torque 
fluctuation control region in this way, the amount of fluctuation of the torque must be detected and the torque must be 
detected to detect the amount of fluctuation of the torque. 

Various methods have been proposed in the past for calculating tiie output torque of each cylinder. That is. if it were 

4(7 possitsle to calculate the output torque of each cylinder, the thus calculated output torque could be used not only to ena- 
ble the control of the lean air-fuel ratio as explained at)Gve. but also other various types of control, so funding the best 
method for calculating the output torque of each cylinder has great significance. Therefore, in the past, various methods 
have been proposed for calculating the output torque of each cylinder. Merrtionir^g typical ones, there are tiie method of 
providing a combustion pressure sensor in the combustion chamber and calculating the output torque based on the out- 

45 put signal of this combustion pressure sensor or the method of, as explained at the beginning, calculating the output 
torque from the difference t^etween tiie square of the first angular velodty <oa and the square of the second angular 
velocity Qob. 

If a combustion pressure sensor is used, there is the advantage tiiat the torque generated in the cylinder provided 
with the combustion pressure sensor can be reliably detected, but there is conversely tiie disadvantage that the com- 

50 txjstion pressure sensor is required. As opposed to this, the angular velocities coa and cob can be calculated from tiie 
output signal of the crank angle sensor provided in the intemal combustion engine in the past, so when calculating the 
output torque based on the angular velocities a>a and cob. there is the advantage that there is no need to mount a new 
sensor. However, in this case, as explained at the beginning, if the engine drive system experiences torsional vibration, 
there is the problem that the generated torque can no longer be accurately detected. It is clear, however, that if this 

55 problem could just be solved, then the method of calculation of the torque based on the angular velocity, which does 
not require a new sensor, would be preferable. Therefore, tiie present invention calculates the generated torque based 
on the angular velocity and thereby can accurately detect the generated torque even if the engine drive system experi- 
ences torsional vibration. 

Next, the new method according to the present invention for calculating tiie drive force generated by each cylinder 
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and the torque generated In each cylinder will be explained. 

Rrst the method of calculating the drive force generated by each cylinder and the torque generated in each cylin- 
der will be explained while referring to Rgs. 6A and 6B showing the steady operating state where the engine drive sys- 
tem does not experience torsional vibration. As explained above* the crank angle sensor 14 produces an output pulse 

5 each time the crankshaft rotates by 30'' crank angle. Further, the crank angle sensor 1 4 is arranged to generate an out- 
put pulse at the top dead center TDC of the compression stroke of the cylinders #1 . #2, #3, and #4. Therefore, the crank 
angle sensor 14 produces an output pulse for each 30"* crank angle from the top dead center TDC of the compression 
stroke of the cylinders #1. #2. #3, and #4. Note that, the ignition sequence of tfie internal combustion engine used in 
the present inventk>n is 1-3-4-2. 

10 In Rgs. 6A and 6B, the vertical axis T30 shows the elapsed time of 30^ crank angle from when the crank angle sen- 
sor 14 produces an output pulse to when it produces the next output pulse. Further. Ta(i) shows the elapsed time from 
the top dead center of the compression stroke (hereinafter refenred to as TDC) to 30** after top dead center of the conv 
pression stroke (hereinafter referred to as ATDC) of the No. i cylinder, while Tb(i) shows the elapsed time from ATDC 
60"* to ATDC 90** of the No. i cylirxier. Therefore, for exanrple. since Ta(1) shows the elapsed time from TDC to ATDC 

IS 30<' of the No. 1 cylinder, while Tb(1) shows the elapsed time from ATDC 60"* to ATDC 90"* of the No. 1 cylinder. On the 
other hand, if the 30^^ crank angle is divided by the elapsed time T30, the result of the division shows the angular velocity 
Qx In this emt)odiment according to the present invention. 30° crank angle/Ta(0 is referred to as the first angular velocity 
ccta in the No. i cylinder, while 30** crank ang!e/TB(i) is referred to as the second angular velocity cob in the No. i cylinder. 
Therefore. 30* crank angle/Ta (1) shows the first angular velocity <oa of the Na 1 cylinder, while 30" crank angle/Tb(1) 

20 shows the second angular velocity cub of the No. 1 cylinder. 

Taking note of the No. 1 cylinder of Rgs; 6A and 6B, when the combustion is started and the combustion pressure 
rises* the elapsed time ^Is from Ta(1) to Tb(1) then rises from Tb(1 ). In other words, the angular velocity of the crank- 
shaft Q> rises from the first angular velocity coa to the second angular velocity cob; then falls again from the second angu- 
lar velocity cGb. That is. the combustion pressure causes the angular velocity oe> of the crankshaft to increase from the 

25 first angular velocity coa to the second angular velocity cob. Rg. 6A shows the case where the comtxistion pressure is 
relatively high, while Fig. 68 shows the case where the combustion pressure is relatively low. From Figs. 6A and 68. 
when the combustion pressure is high, the amount of reduction of the elapsed time (TaO)'Tb(i)) becomes larger com- 
pared with when the conrinistion pressure is low. therefore the amount of increase (cob-coa) of the angular velocity <o 
becomes larger. When the comlxistion pressure becomes higher, the drive force generated by the cylinder becomes - 

30 larger, therefore if the amount of increase (cob-coa) of the angular velocity cd becomes larger, the drive force generated 
by the cylinder becomes larger. Therefore, it is possible to calculate the drive force generated by a cylinder from the dif- ^ 
ference (oob-<na) between the first angular velocity coa and the second angular velocity cob. 

On the other hand, if the moment of inertia of rotation of the engine is 1. the comfcxjstion pressure causes the kinetic 
energy to increase from (1/2) Icoa^ to (1/2) Icob^. The amount of increase of the kinetic energy (1^) *. I ' (cob^-ooa^) 

35 expresses the torque generated by that cylinder, therefore it becomes possible to calculate the torque generated by a ^ 
cylinder from the difference (cob^-coa^) between the square of the first angular velocity coa and the square of the second 
angular velocity cob. 

By detecting the first angular velocity coa and the second angular velocity cob in this way. it becomes possible to cal- 
culate the drive force generated by the corresponding cylinder and the torque generated by the corresponding cylinder 

40 from the detection values. Note that the change in the elapsed time T30 shown in Figs. 6A and 6B differs somewhat 
depending on the engine, therefore the crank angle range for detecting the first angular velocity coa and the crank angle 
range for detecting the second angular velocity cob are set in accordance with the engine so that (cob-coa) best 
expresses the drive force generated t^y the engine or so that (cob^-coa^) best expresses the torque generated by the 
engine. Therefore, depending on the engine, the crank angle range for detecting the first angular velocity coa may be 

45 from before the top dead center of the compression stroke BTDC 30** to TDC. while the crank angle range for detecting 
the second angular velocity cob may be from ATDC 90** to ATDC 120"*. 

Therefore, explaining the method of detecting the angular velocities coa and cob in general terms, the first crank 
angle range is set in the crank angle region from the end of the compression stroke to the beginning of the expansion 
stroke, the second crank angle range is set in a crank angle region in the middle of the expansion stroke a predeter- 

50 mined crank angle away from the first aank angle range, the first angular velodty coa of the crankshaft in the first crank 
angle range is detected, and the second angular vekx;ity of the crankshaft cob in the second crank angle range is 
detected. 

As explained above, by detecting the angular velocities coa and db, it becomes possible to calculate the drive force 
and the torque generated by a corresponding cylinder based on the detected values. The engine drive system, hcw- 
55 ever, experiences torsional vibration occumng at the natural frequency of the drive system due to the successive explo- 
sive actions of the cylinders. If the engine drive system experiences such torsional vibration, it no longer becomes 
possible to accurately detect the drive force and the torque generated by a cylinder t)ased on the angular velocities coa 
and cob. Next, this will be explained with reference to Fig. 7 and Fig. 8. 

Rg. 7 shows the changes in the elapsed time Ta(i) successively calculated for each cylinder when the engine drive 
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system eocperiences torsional vibration. When the engine drive system experiences torsional vibration, this torsional 
vibration causes the angular velocity of a crankshaft to be cyclically increased and decreased, so the elapsed time Ta(i) 
increases and decreases cyclically as shown in Rg. 7. 

On the other hand, Fig. 8 shows the portion where the elapsed time Ta(i) is reduced in an ^larged manner. As 

5 shown in Rg. 8, the elapsed time Ta(l) falls by ho between Ta(1) and Ta(3). This reduction of ho is believed to be due to 
an increase in the amount of torsion due to the torsional vibration. In this case, the amount of decrease of the elapsed 
time due to the torsional vibration between Ta(1) and Ta(3) is believed to increase substantially linearly along with the 
elapse of time, therefore this amount of decrease of the elapsed time due to the torsional vibration is shown by the dif- 
ference between the broken line connecting Ta(1) and Ta(3) and the horizontal line passing through Ta(1). Therefore, 

10 between Ta(1) and Tb(1), the torsional vibration causes the elapsed time to fall by exactiy h. 

That is. Tb(1) is lower in elapsed time than Ta(1). but this lower elapsed time includes the amount of decrease f of 
the elapsed time due to the oonik)ustion pressure arvj the amount of decrease h of the elapsed time due to the torsional 
vibration. Therefore, to find just the elapsed time Tb*(1 ) decreased due to the combustion pressure, it becomes neces- 
sary to add h to Tb(1). That is, when the elapsed time Ta(i) decreases between cylinders (Ta(1)-»Ta(3)). to find just the 

IS elapsed time Tb'(l) decreased due to the combustion pressure, the detected elapsed time Tb(1) ^ corrected in 
the upward direction. In other words, when tiie first angular velocity coa increases between cylinders, the second angu- 
lar velocity cob of the cylinder where the combustion was first performed must be corrected In the downward direction. 

As opposed to this, when Ta(3) increases witii respect to Ta(2), the elapsed time Tb(1) reduced from Ta(1 ) includes 
the amount of decrease of the elapsed time due to the comtxistion pressure and the eimount of increase of tfie elapsed 

20 time due to the torsional vibration. Therefore, in this case, to find just the elapsed time Tb*(1) reduced due to tiie com- 
bustion pressure, the amount of increase of the elapsed time due to the torsional vibration must be subtracted from 
Tb(1). That is, when the elapsed time Ta(i) increases between cylirKlers. to find just the elapsed time Tb'(1) decreased 
due to the combustion pressure, the detected elapsed time Tb(l ) must be corrected in the downward direction. In other 
words, when the first angular velocity coa decreases between cylinders, the second angular velocity cob of tine cylinder 

25 where the combustion was first performed must be corrected in the upward direction. 

As explained above, correction of the second angular velocity db enables the drive force generated t>y each cylin- 
der to be accurately detected from the difference {db-wa) between the first angular velocity <ea ajnd the secorxl angular 
velocity cob even when the engine drive system experiences torsional vibration and enables the torque generated by 
each cylinder to be accurately calculated from the difference (cob^-coa^ between the square of the first angular velocity 

30 coa and the square of the second angular velocity cob. Note, however, that if there is variation in the spaces between the 
outer teeth formed along the outer periphery of the rotor 1 3 (Rg. 1 ), even if the second angular velocity cob is corrected 
as explained above, the drive force and the torque generated by a cylirxler cannot be accurately detected. Next, this will 
be explained with reference to Rg. 9. 

Rg. 9 shows the case where tiie space between the outer tooth of the rotor 13 showing the TDC of the No. 1 cyl- 

35 inder #1 and the outer tooth of the rotor 31 showing ATDC 30" is smaller than the space between other outer teeth. In 
this case, as will be understood from a comparison of Fig. 8 and Fig. 9, the elapsed time Ta(1) will end up becoming 
smaller than the correct elapsed time for 30'' crank angle. Further, at this time, as will be understood from a comparison 
of Rg. 8 arxi Rg. 9, the amount of decrease h' of the elapsed time due to the torsional vibration will end up becoming 
smaller than the correct amount of decrease h and therefore the value off Tb'(1) expressing just the elapsed time 

40 decreased due to the combustion pressure will also erKi up becoming smaller than the correct value. 

Therefore, in this emkxxliment according to the present invention, the ratio KTafi) (= Ta(i)m/Ta(i)) between the 
mean value Ta(i)m of the elapsed times Ta(i) of all of the cylinders and the elapsed time Ta(i) of each cylinder and the 
ratio KTb(0 (= Tb(i) m/Tb(i)) between the mean value Tb(i)m of the elapsed times Tb(i) of all of the cylinders and the 
elapsed time Tb(i) of each cylinder at the time of when the supply of fuel has been stopped in a deceleration operation 

45 when the engine drive system does not experience torsional vibration are found. When the fuel is being supplied, the 
actually detected elapsed time Ta(i) for each cylinder is multiplied by the ratio KTa(i) so as to find the final elapsed time 
Ta(i) for each cylinder and the actually detected elapsed time Tb(i) for each cylinder is multiplied by tiie ratio KTb(i) so 
as to find the final elapsed time Tb(i) for each cylinder. 

Therefore, for example, when, as explained atxjve, the elapsed time Ta(1) actually detected for the No. 1 cylinder 

so #1 is shorter ttian tiie correct elapsed time, the ratio KTa(1) becomes larger than even 1.0, therefore the final elapsed 
time Ta(1) obtained by multiplying tiie actually detected elapsed time Ta(1) with tiie ratio KTa(1) becomes considerably 
close to the correct elapsed time Ta(l). Further, by finding the amount of decrease of the elapsed time due to tiie tor- 
sional vibration h based on the thus obtained final elapsed time Ta(i), this amount of decrease h substantially matches 
with the correct anrount of decrease, therefore even the value of Tb*(1) expressing just the elapsed time due to the com- 

55 bustion pressure substantially shows the cx)nrect value. In this way. in this embodiment according to the present inven- 
tion, even if there is a variation in the spaces between the outer teeth of the rotor 13, the drive force and tiie torque 
generated at each cylinder can be accurately detected. 

On tiie otiier hand, the Ta(i) for each cylinder fluctuates when the vehicle is traveling over a bumpy road. Furtiier, 
at tills time, sometimes the amount of fluctuation of Ta(i) becomes extremely large. Rg. 10 shows the fluctuation in Ta(i) 
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when the vehicle is traveling over a bumpy road. AMP of Rg. 10 shows the difference between the minimum Ta(i) and 
maximum Ta(i), that is. the amplitude. When this amplitude AMP is small, by calculating h shown in Rg. 8 by the method 
explained up to now, the value of Tb*(i) expressing only the elapsed time due to the combustion pressure can be accu- 
rately detected. 

5 However, when the amplitude AMP becomes large, the drive force or the torque generated at a cylinder at which 

Ta(i) becomes maximum or minimum can no longer be accurately detected. That is. in Fig. 1 0. when for example the 
cylinder giving the maximum Ta(i) is first the No. 1 cylinder, the amount of decrease h due to the torsional vibration for 
calculating the Tb'(1) of the No. 1 cylinder #1 is found from the inclination of the broken line connecting Ta(1) and Ta(3) 
in Rg. 10. However, near when the No. 1 cylinder #1 reaches TDC, the amount of increase or the amount of decrease 

10 of the elapsed time due to the torsional vibration changes k>y the smooth curve passing through Ta(2). Ta(1). and Ta(3). 
Therefore, if the value of the amount of decrease h of the No. 1 cylinder #1 with respect to Tb(1) is found from the incli- 
nation of tiie broken line connecting Ta(1) and Ta(3). the value of tiiis amount of decrease h is calculated considerably 
larger than even the actual value. As a result Tb*(1) no longer shows the correct value, therefore the drive force and the 
torque generated at the cylinder can no longer be accurately detected. When ttie amplitude AMP becomes large, tiie 

15 same thing occurs at the cylinder giving the mininmjm Ta(l)* 

Further, in a cylinder where TaO) changes sharply from the Tap) of the cylinder at which the combustion was per- 
formed one time before, the value of h deviates from the actual value, therefore the drive force arKJ tiie torque generated 
by the cylinder can no longer be accurately detected. Therefore, in this embodiment according to the present invention, 
when the amplitude AMP Is large, tiie drive force or the torque for the cylinder at which Ta(l) becomes maximum or mln- 

20 imum is not sought. Further, the drive force or the torque is also not sought for a cylinder where the Ta(i) rapidly changes 
from the Ta(i) for the cylirxier at which the conribustion was perfbnried one ti 

Next, the routines for finding the torque generated at each cylinder will be explained refem'ng to Rg. 1 1 to Rg. 21 . 
Note that, Rg. 21 shows the timing for calculation of the various values performed in each routine. 

Rg. 1 1 shows an interruption routine performed at every SO"" crank angle. Referring to Rg. 1 1 . first the routine (step 

25 100) for calculating tiie elapsed times Ta(i) arKi Tb(i) Is proceeded to. This routine is shown in Rg. 12. Next, the routine 
(step 200) for checking If calculation of the torque Is permitted or not Is proceeded ta This routine is. shown in Rg. 13 
to Rg. 15. Next, the routine for calculating the torque (step 300) Is proceeded to. This routine is shown in Fig. 17, Next, 
the routine for calculating tiie ratios KTa(i) and KTb(i) (step 400) is proceeded to. This routine is shown In Rg. 18 and 
Rg. 19. Next, the routine for processing of tiie counter CDLNIX used for calculation of the torque fluctuation value is 

30 proceeded to. This routine is shown in Rg. 20. 

Referring to Fig. 12 showing the routine for calculation of the elapsed times Ta(i) and Tb{\), first at step 101. the 
time is made the TIMEO. ITie electronic control unit 20 is provided with a free run counter showing the time. The time : 
is calculated from the count value of this free run counter. Next, at step 102. the current time is fetched. Therefore, the 
TIMEO of step 101 expresses the time of 30*> crank angle before. 

35 Next, at st^ 103. whettier the Na 1 cylinder is currently at ATDC 30'' or not is judged. When the No. 1 cylinder is 
not currentiy at ATDC 30**, step 106 is junped to. where whether the No. 1 cylinder is currentiy at ATDC 90*" or not is 
judged. When the No. 1 cylinder Is not currentiy at ATDC 90*", the routine for calculation of the elapsed times Ta(i) and 
Tb(i) is ended. 

As opposed to this, when it is judged at step 103 tiiat tiie Na 1 cylinder is currentiy at ATDC 30**, step 104 is pro- 
40 ceeded to, where the final elapsed time Ta(i) from TDC to ATDC 30*" of the Na 1 cylinder is calculated based on the 
following equation: 

Ta(i) = KTa(i) • (TIME-TIMEO) 

45 That Is, for example, when tiie Na 1 cylinder #1 is currentiy at ATDC 30^ the final elapsed time Ta(l ) from TDC to ATDC 
30* of tiie Na 1 cylinder #1 is calculated from KTa(1) • (TIME-TIMEO). Here. (TIME-TIMEO) expresses the elapsed time 
Ta(1) actually measured from tiie crank angle sensor 14 and KTa(1) is a ratio for conrection of the error due to the 
spaces of the outer teeth of tiie rotor 13. therefore the final elapsed time Ta(l) otstained by multiplying (TIME-TIMEO) 
with KTa(1) comes to accurately express tiie elapsed time when the crankshaft rotates by 30** crank angle. 

so Next at step 105, the flag XCAL(i-1) of the No. (i-1) cylinder where combustion had been performed one time 
before showing that ttie generated torque should be calculated Is set (XCAL(i-1) "1"). In this embodiment according 
to the present invention, as explained above, since the ignition sequence is 1-3-4-2, when the No. 1 cylinder #1 is cur- 
rently at ATDC 30**. the flag XCAL(2) of tiie No. 2 cylinder #2 where tiie combustion had been performed one time 
before shewing that the generated torque should be calculated is set In the same way. when tiie final elapsed time 

55 Ta(3) is calculated as shown in Rg. 21 . tiie flag XCAL (1 ) is set. when the final elapsed time Ta(4) Is calculated, tiie flag 
XCAL (3) is set and when the final elapsed time Ta(2) is calculated, the flag XCAL (4) is set. 

On the other hand, when it is judged at step 1 06 that tiie No. 1 cylinder is currentiy at ATDC 90**. step 1 07 Is pro- 
ceeded to, where tiie final elapsed time Tb(i) from ATDC eO"" to ATDC 90'' of the No. 1 cylinder is calculated based on 
the following equation: 
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TbCO = KTbp) • (TIME-TIMEO) 

That is, for example, when the No. 1 cylinder #1 is cunrently at ATDC 90", the final elapsed time Tb(1) from ATDC 
60*» to ATDC 90** of the No. 1 cylinder #1 is calculated from KTb(1) • (TIME-TIMEO) , In this case as well, since the ratio 
5 KTb(1) for correcting the en-or due to the spaces of the outer teeth of the rotor 13 is multiplied with (TIME-TIMEO), the 
final elapsed time Tb(1) accurately expressed the elapsed time in the period when the crankshaft rotates by 30"* crank 
angle. 

Next, thQ routine for checking permission for calculation of the torque shown in Rg. 13 to Rg. 15 will be explained 
referring to Rg. 16. This routine is set so as to prohibit the calculation of the torque for a specific cylinder when the 
10 amplitude AMP of the fluctuation of Ta(0 (Rg. 1 0) becomes large due to the vehicle traveling over a bumpy road. 

That is, referring to Rg. 13 to Rg. 15. first at step 201 , whether one of the cylinders is cunrently at ATDC 30** or not 
is judged. When none of the cylinders is currently at ATDC 30"*, the processing cycle is ended, while when one of the 
cylinders is at ATDC 30**, step 202 is proceeded to. 

At step 202 to step 204. the maximum elapsed time TSOmax when the elapsed time Ta(i) increases and then 
IS decreases is calculated. That is, at step 202, whether the Ta(i) calculated at the routine shown in Rg. 12 is larger than 
the maximum elapsed time TSOmax or not is judged. When T30max > Ta{i), step 205 is jumped to. while when TSOmax 
^ TaO)> st^ 203 is proceeded to, where Ta(i) is made TSOmax. Next, at step 204, the increase flag XMXREC showing 
that Ta(i) is increasing is set (XMXREC <- "1 ^. then step 205 is proceeded to. 

A step 205 to step 207. the minimum elapsed time T30min when the elapsed time Ta(i) decreases and then 
20 increases is calculated. That is, at step 205, whether the Ta(i) calculated by the routine shown in Rg. 12 is smaller than 
the calculated minimum elapsed time T30min or not is judged. When TSOmin < Ta(i), stq^ 208 is jumped to. while when 
T30min ^ Ta{i), step 206 is proceeded to. where Ta(i) is made TSOmin. Next, at step 207, the decrease flag XMNREC 
showing that Ta(i) has decreased is set (XMNREC "1 "), then step 208 is proceeded to. 

At step 208 to step 214, when the amplitude AMP of the fluctuation of Ta(i) (Rg. 10) has exceeded the setting Aq, 
25 the prohibition flag prohibiting the catoulafion of the torque for the cylinder where Ta(i) k>ecomes maximum is set. That 
is, at step 208, whether T30max > Tap) and XMXREC = "1 " or not is judged. When T30max ^ Tafi) or the Increase flag 
XMXREC is reset (XMXREC « "0"), stisp 215 is jumped to, while when TSOmax > Ta(i) and XMXREC = "1", step 209 is 
proceeded to. 

That is, as shown inRg. 16, assume that at the time t^, the elapsed time Ta(1) of the No. 1 cylinder #1 has become 
30 maximum. In this case, in the intenruption routine performed at the time ti in Rg. 1 6, step 202 to step 203 is proceeded 
to. where the Ta(1) is made TSOmax. then; at step 204. the increase flag XMXREC is set. On the other hand, in the inter- 
ruption routine panned at the time t2 of Rg. 16. step 202 to step 205 is jumped to. At this time, since it is judged at 
step 208 that TSOmax > Ta(3) and XMXREC = "1", step 209 is proceeded to That is. step 209 is proceeded to at the 
time t2 when the elapsed time Ta(Q starts to decrease. At step 209. the maximum elapsed time TSOmax is made 
35 TMXREC. Next, at step 210. the maximum elapsed time TMXREC is sutjtracted by the minimum elapsed time TMN- 
REC (found at the later explained step 2 1 6] so as to calculate the amplitude AM P of the fluctuation of Ta(i). Next, at step 
21 1 . the initial value of the minimum elapsed time TSOmin is made Ta(i). Next at step 212. the increase flag XMXREC 
is reset (XMXREC ''0'*). Next, at step 213, whether the amplitude AMP is larger than the setting Aq or not is judged. 
When AMP < Aq. step 21 5 is jumped to. As opposed to this, when AMP s Aq. slep 21 4 is proceeded to, where the torque 
40 calculation prohibition flag XNOCAL is set (XNOCAL <- "1 "). That is. in the interruption routine performed at the time X2 
of Rg. 16. as explained atxjve. the generated torque of the No. 1 cylinder #1 is calculated. Therefore, in this interruption 
routine, when AMP ^ Aq and the torque calculation prohibition flag XNOCAL is set. the calculation of the generated 
torque of the No. 1 cylirider #1 . that is, the calculation of the generated torque of the cylinder giving the maximum Ta(i), 
is prohibited. 

45 At step 215 to step 221 . when the amplitude AMP of the fluctuation of Ta(i) has exceeded the setting Aq. the prohi- 
bition flag prohibiting the calculation of the torque of One cylinder giving the minimum Ta(i) is set TTiat is. at step 215, 
whetiier TSOmin < Ta{i) and XMNREC = "1" or not is judged, when TSOmin ^ Ta(i) or the decrease flag XMNREC is reset 
(XMNREC =s -0"). step 222 is jumped to. while when TSOmin < TaO) and XMNREC = "r. step 216 is proceeded to. 
That is, as shown in Rg. 16. at the time X3, the elapsed time Ta(1) of the No. 1 cylinder #1 Is considered to have 

so become the minimum. In this case, at the interruption routine performed at the time 13. step 205 to step 206 is pro- 
ceeded to. where the Ta(1) is made TSOmin. then at step 207. the decrease flag XMNREC is set. On the other hand, 
the interruption routine perfomied at tiie time t4 of Rg. 1 6 jumps from step 205 to step 208. At this time, since it is judged 
at step 21 5 that TSOmin < Ta(S) and XMNREC = "r , step 21 6 is proceeded to. That is, step 21 6 is proceeded to at the 
time t4 where the elapsed time Ta(i) starts to be increased. At step 216, the minimum elapsed time TSOmin is made 

55 TMNREC. Next, at step 217, the minimum elapsed time ™nREC is subtracted from the maximum elapsed time 
TMXREC wheretjy the anplitude AMP of tiie fluctuation of Ta(i) is calculated. Next, at step 218, the initial value of the 
maximum elapsed time TSOmax is made Ta(i). Next at step 219, the decrease flag XMNREC is reset (XMNREC <- "0"). 
Next, at step 220, whettier the amplitude AMP is larger than the setting Aq or not is Judged. When AMP < Aq. step 222 
is jumped to As opposed to this, when AMP ^ Aq. step 221 is proceeded to, where tiie torque calculation prohibition 
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flag XNOCAL is set (XNOCAL <- "1"). That is, in the interruption routine performed at the time t4 of Fig. 16. the gener- 
ated torque of the No. 1 cyfinder #1 is calculated. Therefore, in this interruption routine, when AMP ^ Aq and the torque 
calculation prohibition flag XNOCAL Is set. the calculation of tiie generated torque of the Na 1 cylinder #1 . that Is, tiie 
calculation of the generated torque of the cylinder where Ta(i) becomes smallest is prohibited. 

5 At step 222 and step 223. the calculation of the torque of a cylinder where tiie elapsed time Ta(i) changes sharply 

Is prohibited. That is. at step 222. whetiier |Ta 0-2)-Ta G- 1 )l is larger than K q • Fa (i-1 )-Ta (i)| or not is judged. Here, tiie 
constant Kq is a value of about 3.0 to 4.0. When it is judged at step 222 that |Ta0'2)- Ta(i-1)| < Kq - [Ta(M)-Ta(i)| . the 
processing routine Is ended, while when it is judged that |Ta(l-2)-Ta(l-1)| ^ Kq • |Ta(i-1)-Ta(i)| . step 223 is proceeded 
to. where the torque calculation prohibition flag XNOCAL is set. 

10 That is. when tiie interruption routine is at the time ta of Fig. 16. whether at this time rTa(4)-Ta(2)| is 
Kq • |Ta(2)-Ta(1)| or not is judged. As shown in Fig. 16, when Ta(2) changes rapidly from Ta(4). |Ta(4)-Ta{2)| becomes 
larger than Kq * |Ta(2)-Ta(1)| . At this time, the torque calculation prohibition flag Is set and the calculation of the torque 
of the Na 2 cylinder #2 where the elapsed time Ta(i) has changed sharply is prohibited. 

Next, the torque calculation routine shown in Rg. 17 will be explained. Referring to Fig. 17. first, at step 301. 

IS whether the flag XCAL(i-1) showing that the generated torque of tiie No. (i-1) cylinder where contxistion had been per- 
formed one time before should be calculated is set or not is judged. When the flag XCAL(i-1) = "0", that is when the flag 
XCAL(i-1) is not set. ttie processing cycle is ended. As opposed to this, when the flag XCAL(i-1) = "I", that is. the flag 
XCAL(i-l) is set. step 302 is proceeded to. where the flag XCAL(i-1) is reset, then step 303 is proceeded to. 

At step 303. whether the prohibition flag XNOCAL prohibiting the calculation of tiie torque for the cylinder at which 

20 the combustion was performed one time k>efore is reset (XNOCAL = "0") or not is judged. When thte prohibition flag is 
set (XNOCAL = "1"). step 310 is proceeded to, where the prohibition flag XNOCAL is reset As opposed to this, when 
prohitMtion flag is reset, step 304 is proceeded ta That is. only when the flag XCAL is set and prohibition flag XNOCAL 
is reset is step 304 proceeded to. 

At step 304. the amount of change h of the elapsed time due to tiie torsional vibration of the engine drive system 

25 (Rg. 8) is calculated based on the following equation: 

h = {Ta(i-1)-Ta(i)} • 60/180 

That is; as will be understood from Rg. 8. the amount of change h of the elapsed time becomes one-third of' 
30 h^ (= Ta(i-1)-Ta(i)) . Next at step 305. Tb'(irl) expressing only the elapsed time decreased due to the comtHJStion 
pressure is calculated based on .tiie following equation: . 

,) r TO(i-l) = 7bCi-l)+h 

35 That is, when finding Tb'(1) for tiie No. 1 cylinder #1. h {Ta(1)-Ta(3)} • 60/180 and Tb'(1) = Tb(1)+h . Furtiier. 
when finding lb\3) for the No. 3 cylinder #3. h {Ta(3)-Ta(4)} - 60/180 and TbX3) = Tb(3)+h . Next, at step 306, tiie 
generated torque DN(I-I) of the cylinder at which the comtxistion was performed one time before is calculated based 
on the following equation: 

40 DN(i-l) = a>b^-Q>a^ = (30*/Tb'(i-1))^-(30-/TaCi-1)^ 

This generated torque DN(i-1) expresses the torque after elimination of the effect due to the torsional vibration of 
the engine drive system and the effect due to tiie variation in spaces of the outer teeth of the rotor 13. therefore this 
generated torque DN(i-l) expresses the true torque generated due to the combustion pressure. 
45 Note tiiat when finding tiie drive force GN(i-1) generated by each cylinder, tiiis drive force GN(i-1) may be calcu- 
lated based on the following equation: 

GN(i-1) = (30«/rb'(i-1))-(30'*/Ta(i-1)) 

so When the generated torque DN(i-1) is calculated at step 306, step 307 is proceeded to, where tiie amount of fluc- 
tuation of the torque DLN(i-l) in a single cyde of the same cylinder is calculated teased on the following equation: 

DLN(i-l) = DN(i-1)j-DN(i-1) 

55 Here. DN(i-1 )j expresses the generated torque of tiie same cylinder one cyde (720* crank angle) before for DNp-l ). 
Next, at step 308. whether the amount of fluctuation of the torque DLN(i-l) is positive or not is judged- When 
DLN(i-l) ^ 0 . step 310 is jumped to. where tiie cumulative addition request flag XCDLN(i-l) showing that the amount 
of fluctuation of the torque DLN(i-l) of tiie cylinder at which the combustion was performed one time before should be 
cumulatively added is set (XCDLN(i-1 ) ^ "1 As opposed to tills. When DLN(i-l) < 0 , step 309 is proceeded to, where 
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DLN(i-l) is made zero. Next, at step 310 is proceeded to. Note that, the torque of each cylinder repeatedly rises and 
falls, so in this case to find the amount of fluctuation of the torque, rt is sufficient to cumulatively add either the amount 
of increase of the torque or the amount of decrease of the torque. In the routine shown in Rg. 17, just the amount of 
decrease of the torque is cumulatively added, therefore, as explained above, when DLN(i-l) < 0 . DLN(I-I) is made 
5 zero. 

Next, the routine for calculating the ratios KTa(i) and KTb(0 will be explained referring to Fig. 18 and Rg. 19. 

Referring to Rg. 18 and Rg. 19. first, at step 401 . whether the supply of fuel has been stopped during the deceler- 
ation operation or not. that is. whether the fuel has been cut or not. is judged. When the fuel has not been cut. step 415 
is proceeded to. where the cumulative values LTa(i) and LTb(0 of the elapsed times TaO) and Tb(i) are cleared, then the 
10 processing cycle is completed. As opposed to this, when the fuel has been cut. step 402 is proceeded to. where 
whether the amplitude AMP calculated in the routine for checking permission for calculation of the torque is larger than 
the setting Bq or not is judged. When AMP > Bq, step 41 5 is proceeded to. while when AMP ^ Bq. step 403 is proceeded 
to. 

At step 403 to Step 408. KTa(i) is calculated. That is, at step 403. the corresponding elapsed time TaO) each cyl- 
IS inder Is added to the cumulative value LTaO). For example, Ta(1) is added to £Ta(1) and Ta(2) is added to £Ta(2). Next, 
at step 404. whether the Ta(i) for each cylinder has been cumulatively added n number of times each or not is Judged. 
When not cumulatively added n number of times each, step 409 is jumped to, while when cumulatively added n number 
of times, step 405 is proceeded to. At step 405. the mean value Ma (= {i:Ta(1)+i:Ta(2>4-X:Ta(3)-iri:Ta(4)}/4) of the 
cumulative values LTa(0 of the cylinders is calculated. Next, at step 406. the correction value a(i) (= Ma/LTa Q)) for the 
so cylinders is calculated. Next, at step 407. the ratio KTaO) is updated based on the following equation: 

KTa(0 KTa(i)+{a(i)-KTa(i)}/4 

In this way. the ratios KTa(1). KTa(2), KTa(3). and KTa(4) for the cylinders are calculated. For example, if a(1) has 
25 become larger than the KTa(1) used up to then, one-quarter of tiie difference between a(1) and KTa(1) {a(1)-KTa(l)} is 
added to KTa(1), therefore KTa(1) gradually approaches a(1). At step 407. the KTa(i) for each cylinder is calculated, 
then step 408 is proceeded to. where the cumulative value £Ta(i) for each cylinder is cleared. 

On the other hand, at step 409 to step 414. KTb(i) is calculated. That is. at step 409, the corresponding elapsed 
time IbQ) for each cylinder is added to the cumulative value LTbCQ. For example, Tb(1) is added to ETbCl) and Tb(2) is 
30 added to LTb(2). Next, at Step 410. whether the Tb(i) for each cylinder has each t^en cunrujl^^ ^ 
times or not is judged. When not cumulatively added n number of times each, the processing cydie is ended, while when 
cumulatively added n number of times, step 41 1 is proceeded to. At step 41 1 . the mean value 
Mb (^ {£Tb(1)'fLTb(2)-i>LTb(3)+£Tb(4)}/4) of the cumulative values ETbO) Of tiie cylinders is calculated. Next^atstep 
41 2. the correction value p(i) (a Mb/LTbO)) for each cyiirKler is calculated. Next, at step 41 3, the ratio KTb(i) is updated 
35 tiased on the following equation: 

KTb(0 4- KTb(i)+{p(i)-KTb(i)}/4 

In this way. the ratios KTb(1), KTb(2). KTb(3), and KTb{4) for the cylinders are calculated. For exanple. assuming 
40 that p(1) has become larger tiian the KTb(1) used up to then, one-quarter of the difference t>etween p(1) and KTb(1) 
{p(1)-KTb(1)} is added to KTb(l). therefore KTb(l) gradually approaches P(1). When tiie KTb(i) for each cylinder is cal- 
culated at step 413. step 414 is proceeded to. where the cumulative value 2:Tb(i) for each cylinder is cleared. 

Next, the processing of the counter CDLNIX will be explained referring to Fig. 20. The count value of the counter 
CDLNIX is used for the later explained calculation of the torque fluctuation value. 
45 Referring to Rg. 20. first, whether the No. 3 cylinder #3 is currently at ATDC 30"* or not is judged. When the Na 3 
cylinder #3 is cun^entiy not at ATDC 30'', the processing cycle is ended, while when the No. 3 cylinder #3 is currently at 
ATDC 30**, step 502 is proceeded to. At step 502. whether the conditions for calculating the torque fluctuation value 
stand or not is judged. For example, when the conditfons for making the air-fuel ratio lean do not stand or the amount 
of change per unit time APM of the ak>solute pressure of the surge tank 3 is more than the setting or the amount of 
so change per unit time AN of the engine speed is more than a setting, it is judged that the conditions for calculating tiie 
fluctuation value do not stand, while at other times it is judged that the corxiitions for catoulating the fluctuation value 
stand. 

When it is judged at step 502 tiiat tiie conditions for cafoulating the fluctuation value stand, step 508 is proceeded 
to, where the count value CDLNIX is incremented by exactiy 1 . The increment action of this count value CDLNIX Is per- 
55 formed every time the No. 3 cylinder #3 reaches ATDC 30*. that is. every 720** crank angle. Next, at step 509, the aver- 
age value of the engine speed N/^ve ^® average value PM/^ve ^ absolute pressure in the surge tank 3 in the 
period from when the increment action of the count value CDLNIX is started to when the count value CDLNIX is cleared 
are calculated. 

On tiie other hand, when it is judged at step 502 that the conditions for calculating the fluctuation value do not 
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stand, step 503 is proceeded to, where the count value CDLNIX is cleared. Next, at step 504, the cumulative value DLNI 
(i) of the torque fluctuation value DLN(i) for each cylinder (this cumulative value is calculated by the later explained rou- 
tine) is cleared. Next, at step 505, the cunrulative count value CDLNI (i) for each cylinder (this cumulative count value 
is calculated by the later explained routine) is cleared. 

5 Next, at step 506. the target torque fluctuation value LVLLFB is calculated. In this embodiment according to the 

present invention, as explained later, the air-fuel ratio is feedback controlled so that the calculated torque fluctuation 
value becomes this target torque fluctuation value LVLLFB. This target torque fluctuation value LVLLFB, as shown by 
Rg. 22 showing the equivalent fluctuation value l3y the solid line, becomes larger the higher the absolute pressure PM 
in the surge tank 3 and becomes larger the higher the engine speed N. This target torque fluctuation value LVLLFB is 

10 stored in the ROM 22 in advance in the form of a map shown in Fig. 22B as a function of the absolute pressure PM in 
the surge tank 3 and the engine speed N. Next at step 507, the mean torque fluctuation value DLNlSM(i) of each cyl- 
inder (tiiis mean torque fluctuation value is calculated by the later explained routine) is made the target torque fluctua- 
tion value LVLLFB calculated from the map of Rg. 22B. 

Rg. 24 shows the repeatedly executed main routine. In this main routine, first the routine for calculation of the 

IS torque f hictuation value (step 600] is executed. This routine is shown in Rg. 25 and Fig. 26. Next, the routine for calcu- 
lation of tiie lean limit feedback correction coefficient FLLFB (step 700) is executed. This routine is shown in Rg. 27. 
Next, when the predetermined crank angle is reached, the routine for calculation of the injection time (step 800) ^ exe- 
cuted. This routine is shown in Rg. 28. Next, the other routines (step 900) are executed. 

Next, the routine for calculation of the torque fluctuation value shown in Rg. 25 and Rg. 26 will be explained. 

20 Referring to Fig. 25 and Rg. 26, first, at step 601 , whether ttie cumulative addition request flag XCDLN(i) showing 
that the amount of fluctuation of the torque DLN0) shouki be cumulatively added is set (XCOLN(i) » "1 ") or not is judged. 
When the cumulative addition request flag XCDLN(i) is not set, step 609 is jumped to. while when the cumulative addi- 
tion request flag XCDLN(i) is set, step 602 is proceeded to. At step 602, the cumulative addition request flag XCDLN(i) 
is reset. Next, at step 603, the amount of fluctuation of the torque DLN(i) is added to the cumulative value OLNI(i) of the 

25 amount of fluctuation of the torque. Next, at step 604. the cumulative count value CDLNI(i) is incremented by exactty 1 . 
That is. for example, at step 601 . if tiie cumulative addition request flag XCDLN(I) is set for tiie No. 1 cylinder, this flag 
XCDLN(1) is reset at step 602, the amount of fluctuation of the torque cumulative value DLNI(1) is calculated at step 
603, and the cumulative count value CDLNI(I) is incremented by exactly 1 at step 604. . 

Next, at step 605. whether the curhulative count value CDLNI(i) has become "8" or not is judged: When CDLNI(i) 

30 Is not '*8'*^ step 609 is jumped to, while when CDLNI(i) becomes "8", step 606 is proceeded to^ where the cumulative ' 
value DLNI(i) of the amount of fluctuation of tiie torque is cleared. Next, at step 607, the cumulative count value 
CDLNI(i) is reset. Next, at step 608, the mean torque fluctuation value DLNISM(i) is cafoulated from the following equa- 
tion: 

35 DLNiSM(i} = DLNISM(i)+{DLNI0)-DLNISM(i)}/4 

That is, when there is a difference between the calculated amount of fluctuation of the torque cumulative value 
OLNI(i) and the previously used mean amount of fluctuation of the torque DLNISM(i). the value of the difference 
{DLNI(i)-DLNISM(i)} multiplied by 1/4 is added to the mean aniount of fluctuation of the torque DLNISM(i}. Therefore. 

40 for example^ when the cumulative count value CDLNi(l) for the Na 1 cylinder #1 becomes '*8". at step 606. tiie mean 
torque fluctuation value DLNISM(l) is calculated. 

Next, at step 609, whether the count value CDLNIX calculated at the routine shown in Rg. 20 has become "8** or 
not is judged. When CDLNIX is not "8", the processing cyde is ended, while when CDLNIX becomes "8", st^ 610 is 
proceeded to. where tiie mean value DLNISM (= {DLNISM (1 )+DLNISM (2)+DLNlSM (3)+DLNISM (4)}/4) of the mean 

45 torque fluctuation values DLNISM (i) of the cylinders is calculated. Next, at step 61 1 . the count value CDLNIX is cleared. 
In this way, the value DLNISM expressing the amount of fluctuation of the torque of the engine is calculated. 

Note that as explained above, the count value CDLNIX is incremented by exactiy 1 with each 720*" crank angle. 
Unless the calculatfon of the torque is prohitxted for one of the cylinders, when the count value CDLNIX has become 
"8". the cumulative count values CDLNI(I), CDLNI(2), CDLNI(3). and CDLNI(4) for all of the cylinders have already 

so become "8**. Therefore, in this case, the mean torque fluctuation value DLNiSM(i) for all of the cylinders is calculated. 
On the other hand, for example, if the calculation of the torque for the No. 1 cylinder #1 Is prohibited, when the count 
value CDLNIX has become "8", just the cumulative countvalueCDLNI(l)ofttie No. 1 cylinder #1 does not become '*8^ 
so the new amount of fluctuation of the torque cumulative value DLNI(1) for the No. 1 cylinder #1 is not calculated. 
Therefore, in this case, when finding the mean value DLNISM at step 610, the previously calculated amount of fluctua- 

55 tion of the torque cumulative value DLNISM(I) is used just for the No. 1 cylinder #1 . 
Next, the routine for calculation of FLLFB will be explained referring to Rg. 27. 

Referring to Rg. 27. first, at step 701 , whether the conditions for updating the lean limit feedback correction coeffi- 
dent FLLFB stand or not is judged. For example, at tiie time of engine warmup or when the operating state of the engine 
is not in the learning region enclosed by the broken lines in Rg. 5, it is judged tiiat the conditions for updating do not 
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stand, while at ottier times it is judged that the conditions for updating stand. When the conditions for updating do not 
stand, the processing cyde is ended, while when the conditions for updating stand, step 702 is proceeded to. 

At step 702. tiie target torque fluctuation value LVLLFB is calculated from the absolute pressure PM in tiie surge 
tank 3 and the engine speed N based on the map shown in Fig. 22B. Next at step 703 and step 704. tiie levels of torque 
5 fluctuation LVLH(n) and LVLL(n) shown in the following equatior^ are calculated based on the fluctuation amount judge- 
ment values DH(n) arxJ DL(n) in accordance witii the target torque fluctuation value LVLLFB: 

LVLH(n) = LVLLF&4.DH(n) 

10 LVLUn) =: LVLLFBi.DL(n) 

Here, the fluctuation amount judgement values DH(n) and DL(n) are determined In advance as shown in Rg. 23A. 
That is, as will be understood from Rg. 23A. three positive values are set for DH(n) which are In the relationship of 
DH(3) > DH(2) > DH(1). Furtiier. these DH(1). DH(2). and DH(3) gradually increase as the target toique fluctuation 
IS value LVLLFB becomes larger. On the otiier hand, three negative values are set for DL(n) which are in the relationship 
of DL(1 ) > DL(2) > DL(3). Further, the absolute values of these DL(1). DL(2), and DL(3) gradually increase as the target 
torque fluctuation value LVLLFB becomes larger. 

Assume however that the target toique fluctuation value LVLLFB calculated at step 702 is the value shown by the 
broken line. In this case, at step 703. the values of DH(1). DH(2). and DH(3) on tiie broken tine plus the target torque 
20 fluctuation value LVLLFB are made the levels of torque fluctuation LVLH(1 ). LVLH(2). and LVLH(3) and. at step 704. the 
values of DL(1), DL(2). and DL(3) on the broken line plus the target torque fluctuation value LVLLFB are made the levels 
of torque fluctuation LVLL(1). LVLL(2), and LVLL(3). 

On the other hand, the feecfi^ack correction values +a^, +a2* "^s* ^^i* "^4 ^® determined in 

advarTce for the regions between the levels of torque fluctuation LVLH(n) and LVLL(n) as shown in Rg. 23B. For exam- 
25 pie. the feedt^ck correction value becomes +a2 for the region where the level of torque fluctuation is between LVLH(1 ) 
. and LVLH(2). These feedback con-ection values are +a4 > +a3 > +a2 > +^1 ^ -b^ > -b2 > -ba > -b4. The feedback cor- 
rection values -fa^, +a2> +d3» -bi, - b2. -b^, and -b4 shown in Rg. 23B are shown in the corresponding regions of 
Rg. 23A. 

When the levels of torque fluctuation LVLH(n) and LVLL(n) are calculated at step 703 and step 704, step 705 is pro- 

30. ceeded to, where whether the mean value OLNISM of the torque fluctuation value calculated in the routine for calcula* 
tion of the torque fluctuation value shown in Rg. 25 and Rg. 26 is t>etween the levels of torque fluctuation LVLH(n) and 
LVLL(n) shown in Rg. 23B or not is juJged. Next at step 706. the corresponding feedback correction value DLFB is cal- 
culated. For example, when the target fluctuation level LVLLFB is the value shown by the broken line in Rg. 23A and 
the calculated mean value DUMISM of tiie toique fluctuation value Is ksetween LVLH(1) and LVLH(2) of Rg. 23B, that is 

35 the deviation of tiie mean value DLNISM of the torque fluctuation value from the target fluctuation level LVLLFB is 
between DH(1) and DH(2) on the broken line in Rg. 23A. the feedback correction value DLFB is made +a2- 

Next, at step 707. what lean limit feedback correction coefficient of which learning region shown in Rg. 5 the lean 
limit feedback conrection coeff icient FLLBFij to be updated based on the average value of the engine speed N/^ve 
the average value PM^ve of the absolute pressure in tiie surge tank 3 found at step 509 of the processing routine of 

40 CDLNIX shown in Fig. 20 is determined. Next, at step 708, the lean limit feedback connection coefficient FLLFBij deter- 
mined at step 707 is increased by tiie feedback correction value DLFB. 

That is. as explained above, when for example. DLNISM > LVLLFB and LVLH(1) < DLNISM < LVLH{2), the lean limit 
feedback correction coefficient FLLFBij is increased by +a2. As a result, the air-fuel ratio becomes smaller, so the 
amount of fluctuation of the torque of each cylinder is reduced. On the other hand, when DLNISM < LVLLFB and 

45 LVLL(1) > DLNISM > LVLL(2). the lean limit feedback correction coefficient FLLFBij is increased by -bg. As a result, the 
air-fuel ratio becomes large, so the amount of fluctuation of the torque of the cylinders is increased. In this way the air- 
fuel ratio at the time of lean operation is controlled so that the mean value DLNISM of the amount of fluctuation of tiie 
torque of all of the cylinders becomes tiie target torque fluctuation value LVLLFB. 

Note that, when the conditions for calculation of the torque fluctuation value in the routine shown in Rg. 20 do not 

so stand, at step 507. the DLNISM(i) for all of the cylinders is made LVLLFB and therefore the mean value DLNISM of the 
torque fluctuation value is also made the target torque fluctuation value LVLLFB. Therefore, at this time, the lean limit 
feedback correction coefficient FLLFBij is not updated. 

Next, the routine for calculation of the fuel injection time will be explained with reference to Fig. 28. 

Referring to Rg. 28. first, at step 801 . tiie basic fuel injection time TP is calculated from the map shown in Rg. 2. 

55 Next, at step 802, whether the operating state is one in which a lean operation should be performed or not is judged 
.When the operating state rs one in which a lean operation should be performed, step 803 is proceeded to, where tiie 
value of the stoichiometric air-fuel ratio feedback correction coefficient FAF is fixed at 1 .0, Next, at step 804, the lean 
correction coefficient FLEAN is calculated from the map shown in Fig. 4, then the lean limit feedb>ack correction coeffi- 
cient FLLFB is read from the map shown in Fig. 5. Next, at step 809, tiie fuel injection time TAU is calculated from tiie 
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following equation: 

TAU = TP • FLEAN • FLLFB • FAF+TAUV 

5 As opposed to this, when it is juJged at step 806 that the operating state is not one where a lean operation should 

be perfonmed. that is. when the air-fuel ratio should be made the stoichiometric air-fuel ratio, step 806 is proceeded to. 
where the lean correction coefficient FLEAN is fixed at 1 .0, then, at step 807. the lean limit feedback correction coeff i- 
cient FLLFB is fixed at 1 .0. Next, at step 808. the stoichiometric air-fuel ratio feedback correction coefficient FAF is con- 
trolled based on the output signal of the air-fuel ratio sensor 17 so that the air-fuel ratio becomes the stoichiometric air- 

10 fuel ratia Next, step 809 is proceeded to. where the fuel injection time TAU is calculated. 

By using the method explained at>ove, the drive force generated at each cylinder or the torque generated at each 
cylinder can be accurately detected. Further, when the amount of fluctuation of the drive force or torque is detected from 
the detected drive force or torque, the amount of fluctuation of these drive force or torque can be accurately detected. 
However, a crankshaft experiences natural torsional vibration due to the explosive force successively generated in the 

15 cylinders. Therefore, to accurately detect the drive force generated at each cylinder and the torque generated at each 
cylinder, it is necessary to consider the natural torsional vit>ration of the crankshaft as well. 

Next, an explanation will be made of a second embodiment which enables more accurate detection of the drive 
force generated at each cylinder or torque generated at each cylinder by considering the natural torsional vibration of 
the crankshaft 

20 That is. when an engine is operated, the crankshaft experiences its natural torsional vibration. The amplitude of the 
torsional vibratkui becomes larger the higher the engine speed. This torsional vibration of the crsuikshaft occurs in the 
form of various orders of torsional vibration superposed on each other. Among these, in particular the rotational sixth 
order of torsional vibration (torsional vibration using 60** crank angle as cyde) has a large effect on the detection of the 
first angular velocity caa. Fig. 29 shows the relationship between the amplitude of this rotational sixth order of torsional 
25 vtoration and the engine speed N. As shown in Rg. 29. when tiie engine speed N is low. the anplitude of the rotational 
sixtii order of torsional vibration is small, therefore, at this time, as shown in Rg. 30A. the angular velocity a> rises rela^ 
1 tively smoothly due to tiie combustion pressure when tiie top dead center TDC of tiie compression stroke is exceeded. 

As opposed to this, when the engine speed N becomes high, as shown in Rg. 29. the amplitude of the rotational sixth 
' order of torsional vibration becomes larger. As a result as shown by the arrow mark Z in Rg. 30B, the angular velocity 
30 CD falls sharply for a time due to the effect of tiie rotational sixth order off torsional vibration when the top dead center - 
V TDC of the compression stroke is exceeded; ^ 

Note, however, that if the angular velocity co falls in this way due to tiie rotational sixth order of torsional vibration. 
- the latter half of the crank angle range for detecting the first angular velocity oa will end up overlapping the first half of 
the region of reduction of the angular velocity a> due to the rotational sixUi order of torsional vibration. As a result due 
35 to the effect of this rotational sixth order of torsional vibration, the first angular velocity coa can no longer be accurately 
detected. If the first angular velocity coa can no longer be accurately detected, the drive force and the torque of each 
cylinder calculated t>ased on the first angular velocity coa and second angular velocity cob will no longer show the true 
values. 

In this way. when the engine speed N becomes high, the amplitude of the torsional vibration of the crankshaft 

40 t>ecomes larger and it no longer becomes possible to accurately detect the first angular velocity ooa. At this time, how- 
ever, the amplitude of the torsional vibration generated at the crankshaft as explained earlier, differs depending on the 
crankshaft position. That is, if in Rg. 31 , C is a crankshaft, F is a flywheel. R is a rotational body witii a smaller moment 
of inertia of rotation than the flywheel, for examp>le. a rotor, the amplitude of the torsional vibration generated at the 
crankshaft C will differ according to tiie position of the crankshaft as shown in Rg. 31 . The amplitude of the torsional 

45 vibration of the crankshaft C at the crankshaft position of the No. 4 cylinder #4 closest to tiie flywheel F with a large 
moment of inertia of rotatfon becomes smallest, while the amplitude of the torsional vibration of the crankshaft C at the 
crankshaft position of the No. 1 cylinder #1 furthest away from tiie flywheel F becomes largest. Further, the anrtplitudes 
of torsional vibration generated at the different positions of the crankshaft C become larger the higher tiie engine speed. 
Therefore, as shown in Rg. 29. if the engine speed N becomes higher, tiie crankshaft C experiences torsional 

50 vibration, but at this time, the amplitude of the torsional vibration generated at the crankshaft position of the No. 4 cyl- 
inder #4 is relatively small even if the engine speed is high, therefore it becomes possible to detect the generated torque 
for ttie No. 4 cylinder #4 regardless of the engine speed relatively accurately. However, for example, the amplitude of 
the torsional vibration occurring at the crankshaft position of the Na 1 cylinder #1 becomes considerably large the 
higher tiie engine speed N. therefore for tiie No. 1 cylinder #1, the drive force and the torque generated at the No. 1 

55 cylinder #1 can no longer be accurately detected. 

Therefore, in the second embodiment, tiie calculation of the drive force and the torque is prohibited in the order of 
the cylinders of the crankshaft positions where tiie amplitude of the torsional vibration tsecomes larger as the engine 
speed N becomes higher. Explaining this in more detail, in the second embodiment, when tiie engine speed N is low. 
the drive force and the torque generated at the cylinders are calculated for all of the cylinders #1 , #3, and #4. That 
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is. when the engine speed N is low. the annplitude of the torsional vibration of the crankshaft is small, therefore at this 
time, the torsional vibration of the crankshaft C has almost no effect on the first angular velocity oaa. so the drive force 
and the torque are calculated for all of the cylinders #1 . #2, #3. and #4. 

On the otiier hand, if the engine speed N becomes high, the amplitude of the torsional vibration at tiie position of 

5 the crankshaft at the side away from the flywheel F becomes larger. Therefore, at this time, calculation of the drive force 
and the torque for the No. 1 cylinder #1 and tiie No. 2 cylinder #2 is prohibited. That is. at this time, the drive force and 
the torque are calculated only for the No. 3 cylinder #3 and tiie No. 4 cylinder #4. When the engine speed N becomes 
further higher, the amplitude of the torsional vibration at the crankshaft position near tiie flywheel F becomes larger. 
Therefore, at this time, calculation of the drive force and the torque for the No. 1 cylinder #1 . No. 2 cylinder #2. and No. 

10 3 cylinder itQ is prohibited. That is, at this time, the drive force and the torque for only the Na 4 cylinder #4 is calculated. 
If the calculation of the drive force and the torque for cylinders at crankshaft positions where tiie amplitude of the 
torsional vibration is large is prohibited in this way. the risk of calculating a mistaken drive force and the torque is elim- 
inated. 

In the second embodiment, however, the fuel injection time TAU Is calculated based on the following equation: 

IS 

TAU =« TP • FLEAN • FLLFB • KGTP(i) • FAF+TAUV 

Here, as explained atxive. TP shows the basic fuel injection time, FLEAN shows the lean correction coefficient, 
FLLFB shows tiie le£in limit feedback correction coefficient, FAF shows the stoichiometric air-fuel ratio feedback correc- 
20 tion coefficient, and TAUV shows the invalid Injection time. Note that in the above equation, KGTP(i) shows the intercyl- 
inder correction coefficient. The intercylinder correction coefficient KGnrP(i) is the correction coefficient for.correcting 
the amount of fuel injection for each cylinder so that the times required for the expansion stroke of the cylinders #1 . #2. 
#3. and #4 become equal. 

In the secorxl emfcxxJiment when the air-fuel ratio should be maintained at the stoichiometric air-fuel ratio, the lean 
25 correction coefficient FLEAN, lean limit feedback correction coefficient FLLFB, and intercylinder correction coefficient 
KGTP(i) are fixed at 1 .0 and the stoichiometric air-f uei ratio feedk>ack correction coefficient FAF is controlled based the 
output signal of the air-fuei ratio sensor 1 7. On the other hand, when the air-fuel ratio should be maintained at the target 
lean air-fuel ratio, the stoichiometric air-ftiel ratio feedback correction coefficient FAF is fixed at 1 .0, that is. the feedback 
control based on the output signal of the air-fuel ratio sensor 1 7 is stopped, and the air-fuel ratio is controlled to the tar- 
so. get lean air-fuel ratio t>ased on the lean correction coefficient FLEAN, the lean limit feedback correction coefficient ' 
FLLFB. and the intercylinder correction coefficient KGTP(i). 

Next, the routine for firrding the torque generated at each cylinder will be explained referring to Rg. 32 to Rg. 35. ■ 
Rg. 32 shows the interruption routine performed at every 30'' crank angle. Referring to Rg. 32. first, the routine 
(step 11 00) for calculating the elapsed times Ta(i) and Tb(i) is proceeded to. This routine is shown in Rg. 33. Next the 
35 routine (step 1200) for checking if the calculation of the torque is permitted or not is proceeded to. This routine is shown 
in the previously explained Rg. 13 to Rg. 15. Next, tiie routine for calculating the torque (step 1300) is proceeded to. 
This routine is shown in Fig.34 and Rg. 35. Next, the routine for calculating the ratios KTa{i) and KTb(i) (step 1400) is 
proceeded to. This routine is shown in the previously explained Rg.18 and Rg. 19. Next, tiie routine (step 1500) for 
processing of the counter CDLNIX used for calculation of the torque fluctuation value is proceeded to. This routine Is 
40 shown in the previously explained Rg. 20. 

Referring to Rg. 33 showing the routine for calculating the elapsed times Ta(i) and Tb(i), first, at step 1101. the time 
is made the TIMEO. The electronic control unit 2 is provided with a free run counter showing the time. The time is cal- 
culated from the count value of this free run counter. Next, at step 1 102. the current time is fetched. Therefore, the 
TIMEO at step 1101 comes to express the time of 30" crank angle before. 
45 Next, at step 1 103. whether the No. 1 cylinder is currentiy at ATDC 30"* or not is judged. When the No. 1 cylinder is 
not cun-entiy at ATDC 30*. step 1 106 is jumped to. where whether the Na 1 cylinder is cun-entiy at ATDC 90" or not is 
judged. When the No. 1 cylinder is not currently at ATDC 90"*. step 1108 is junped to. 

As opposed to this, when it is judged at step 1 103 that the Na 1 cylinder is cunrentiy at ATDC 30*". step 1 104 is 
proceeded to, where the final elapsed time Ta(i) from TDC to ATDC 30'' of the No. 1 cylinder is calculated based on tiie 
so following equation: 

Tap) = KTa(i) • (TIME-TIMEO) 

That is. for exanple. when the No. 1 cylinder #1 is currentiy at ATDC 30*'. the final elapsed time Ta(1) from TDC to 
55 ATDC 30" of the Nal cylinder #1 is calculated from KTa(1) -(TIME-TIMEO). Here. (TIME-TIMEO) expresses the 
elapsed time Ta(1) actually measured from the crank angle sensor 14 and KTa(1) is a ratio for correcting the error due 
to the spaces between the outer teeth of tiie rotor 13. therefore the final elapsed time Ta(l) obtained by multiplying 
(TIME-TIMEO) with KTa(1) accurately expresses the elapsed time when tiie crankshaft rotates by 30"* crank angle. 

Next, at step 1105. the flag XCAL(i-1) of the Na (i-l) cylinder where combustion had been performed one time 
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before showing that the generated torque should be calculated is set (XCAL(i-1) 4- "1 In this embodiment according 
to the present invention, as explained above, the ignition sequence is 1 -3-4-2. so when the No. 1 cylinder #1 is currerrtly 
at ATDC 30"*. the flag XCAL(2) of the Na 2 cylinder #2 where combustion had been performed one time before showing 
that the generated torque should be calculated is set. In the same way. when the final elapsed time Ta(3) is calculated 
5 as shown in Rg. 21, the flag XCAL(1) Is set, when the final elapsed time Ta(4) is calculated, the flag XCAL(3) is set, 
and when the final elapsed time Ta(2) is calculated, the flag XCAL(4) is set. 

On tiie other hand, when it is judged at step 11 06 tiiat the No. 1 cylinder is currently at ATDC 90*", step 1 1 07 is pro- 
ceeded to. where the final elapsed time Tb(i) from ATDC 60"* to ATDC 90* of the No. 1 cylinder is calculated based on 
the following equation: 

10 

TbG) = KTb(0 • (JIME-TIMEO) 

That is, for example, when the Na 1 cylinder #1 is cunrentiy at ATDC 90**, the final elapsed time Tb(1) from ATDC 
60** to ATDC 90** of the No. 1 cylinder #1 is calculated from KTb(1) * (TIME-TIMEO). In tiiis case as well, since tiie ratio 
IS KTb(1) fbr confection of the error due to the spaces between the outer teetfi of the rotor 13 is multiplied with (TIME- 
TIMEO), the final elapsed time Tb(1 ) comes to accurately express the elapsed time when the crankshaft rotates by 30" 
crank angle. 

Next, at st^ 1 1 08. whetiier the No. 1 cylinder is cunrentiy at ATDC 21 0* or not is judged. When the No. 1 cylinder 
is currentiy at ATDC 210**, tiie elapsed time T180(i) from ATDC 30** to ATDC 210* of the No. 1 cylinder, ttiat is. tiie 

20 elapsed time T1 80(i) of the expansion stroke of the No. 1 cylinder is calculated. While not explained in detail in the sp^- 
ification of this application, the difference between the elapsed time T180(i) of the expansion stroke of a cylinder and 
the elapsed time T1 80(i) of the next expansion stroke of the cylirxJer is calculated for each cylinder and the intercylinder 
correction coefficient KQTP(i) for each cylinder is calculated so tiiat the difference for each cylinder becomes smaller. 
Note that, tiiis intercylinder correction coeffident KGTP(i) is disclosed in detail in for example Japanese Unexamined 

25 Patent Publication (Kokai) No. 4-370346. 

Next, the torque calculation routine shown in Rg. 34 and Rg. 35 will be explained. Refemng to Rg. 34 and Rg. 35. 
first, at step 1301, whether the flag XCAL(i-1) of the No. (i-1) cylinder where the connbustion was performed one time 
before showing that the generated torque should be cafoulated is set or not is judged. When the flag XCALO-1) = "0", . 
that IS. when the flag XCALO-1) is not set. the processing cycle is ended. As opposed to this, when the flag XCAL(i-1) 

30 s "1". tiiat is. when tiie flag XCAL(i-1) is set, step 1302 is proceeded toj where tiie flag XCAL(i-1) is reset then step i 
1303 is proceeded to. . ^ 

At step 1303. whether the engine speed N is lower than the predetermined set speed NL or not is judged. When N 
< NL, step 307 is proceeded to, where whetiier the prohibition flag XNOCAL-prohit)itihg tiie calculation of tiie torque for , 
the cylinder at which the combustion was performed one time before is reset CXNOCAL » *'0'*) or not is judged. When . 

35 this prohibition flag is set (XNQCAL - "1 "), step 1 3 1 5 is proceeded to, where tiie prohibition flag XNOC AL is reset As 
opposed to this, when the prohibition flag is reset step 1308 is proceeded to. That is. only when the flag XCAL is set 
and prohibition flag XNOCAL is reset is step 1308 proceeded to. 

At step 1308. the amount of fluctuation h of the elapsed time (Rg. 8) due to the torsional vibration of the engine 
drive system is caknjlated Isased on the following equation: 

40 

h = {Ta(i-1)-TaG)) • 60/180 

That is. as will be understood from Rg. 8, the amount of fluctuation h of the elapsed time becomes one-third of 
^0 (= Ta(i-1}-Ta(i}) . Next, at step 1309. 7bXi-1) expressing only the elapsed time due to the comtxjstion pressure is 
45 calculated from the following equation: 

TbXi-l)= Tb(i-1)+h 

That is. when finding Tb*(1) for tiie No. 1 cylinder #1. h = {Ta(1)-Ta(3)} • 60/180 and Tb'(1) = Tb{1)+h . Furtiier, 
50 When finding Tb'(3) for ttie No. 3 cylinder #3. h = {Ta(3)-Ta(4)} • 60/180 and Tb'(3) = Tb(3)+h . Next, at step 310. tiie 
generated torque DN(i-1} of the cylinder at which the comtxjstion was performed one time before is calculated based 
on the following equation: 

DN(i-1) = ©b^-a>a^ = (SO^/TbO-IM^-CSO^/Tap-l))^ 

55 

This generated torque DN(i-l) expresses the torque after elimination of the effect of the torsional vibration of the 
engine drive system and the effect of the variation in spaces between the outer teeth of the rotor 1 3. 

Note tiiat when finding ttie drive force GN(i-1) generated by a cylinder, this drive force GN(i-1) may be calculated 
based on the following equation: 
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GNa-1) = (30*/Tb-(i-1)H30*/ra(i-1)) 

At step 1310, the generated torque DN(i-l) is calculated, then step 131 1 is proceeded to, where the amount of fluc- 
tuation of the torque DLN(i-1) in one cyde of the same cylinder is calculated teased on the follcwing equation: 

5 

DLNG-I) = DN(i-1)j-DN(i-1) 

Here. DN(i-1)j expresses the generated torque of the same cylinder one cyde (720' crank angle) before for DNO-I ). 

Next, at step 1312, whether the amount of fluctuation of the torque DLN(i-l) is positive or not is judged. When 
10 DLN(i-1)^ 0. step 1314 is jumped to. where the cumulative addition request flag XCDLN(i-1) of the cylinder at which the 
combustion was performed one time before showing that the amount of fluctuation of the torque DLJSI(i-l) should be 
cumulatively added is set (XCDLN(i-1) ^ "1 "). As opposed to this, when DLN(i~^) < 0* '^313 is proceeded to, where 
the DLN(i-1) is made zero, then step 1314 is proceeded to. Note that the torque of each cylinder repeatedly rises and 
falls. To find the amount of fluctuation of the torque in this case, it is sufficient to cumulatively add the amount of increase 
15 of the torque or the etmount of decrease of the torque. In the routine shown in Rg. 34 and Rg. 35, only the amount of 
decrease of the torque is cumulatively added, therefore, as explained at^ove, when DLN(i-1) < 0, DLN(i-l) is made zero. 

In this way. when the engine speed N is lower than even the set speed NL, so long as the prohibition flag XNOCAL 
is not set. the generated torque DN(i-1) of the cylinders #1 . #2. #3. and #4 and the amount of fluctuation of the torque 
DLN(i-1) are successively calculated. 
20 On the other hand, at step 1303. when it is judged that N ^ NL. step 1304 is proceeded to, where whether the 
engine speed N is lower tfian the predetermined set speed NH (> NL) or not is judged. When N < NH, step 305 is pro- 
ceeded to. where whether the cylinder (i-1) at which combustion was performed one time t)efore was the Na 3 cylinder 
#3 or the No. 4 cylinder #4 or not is judged. When the cylinder {i-1 ) at which combustion was performed one time before 
was the No. 3 cylinder #3 or the No. 4 cylinder #4. step 1307 is proceeded to, while when the cylinder (i-1) at which com- 
25 bustion was performed one time before was neither the No. 3 cylinder #3 nor the No. 4 cylinder #4. the processing cycle 
is ended. Therefore, when the engine speed N is NL ^ N < NH. so long as the prohibition flag XNOCAL is not set. the 
generated torque DNG-1) and the amount of fluctuation of the torque DLN(i-1) for only the Na 3 cylinder #3 and the No. 
4 cylinder #4 are calculated. 

On the other hand, when it is judged at step 1 304 that N ^ NH. step 306 is proceeded to, where whetho- the cylinder 
- 30 (i-1) at which combustion was performed one time before was the No. 4 cylinder #4 or not is judged. When the cyllrKler 
(i-1) at which combustion was performed one time before was the No. 4 cylinder #4, step 1 307 is proceeded to, while 
when the cylinder (i-1 ) at which comtxistfon was performed one time t>efore was not the No. 4 cylinder #4. the process- 
ing cyde is ended. Therefore, when the engine speed N is N ^ NH, so long as the prohibition flag XNOCAL is not set. 
the generated torque DN(i-l ) and the amount of fluctuation of the torque DLN(i-1 ) for only the No. 4 cylinder #4 are cal- 
35 culated. 

Fig. 36 shows the repeatedly executed main routine. In this main routine, first, the routine for calculation of tiie 
torque fluctuation value (step 1600) is executed. Tliis routine is shown in Fig. 37 and Fig. 38. Next the routine for cal- 
culation of the lean limit feedtiack correction coeff ident FLLFB (step 1 700) is executed .This routine is shown in the pre- 
viously explained Rg. 27. Next, when the predetermined crank angle is reached, the injection time calculation routine 
40 (step 1800) is executed. This routine is shown in Rg. 39. Next, the otiier routines (step 1900) are executed. 

Next, the routine for calculation of the torque fluctuation value will be explained refemng to Rg. 37 and Rg. 38. 
Referring to Fig. 37 and Fig. 38. first at step 1 601 . whether the cumulative addition request flag XCDLN(i) showing 
that the amount of fluctuation of the torque DLN(i) should be cumulatively added is set (XCDLN(i) s "1 ") or not is judged. 
When the cumulative addition request flag XCDLN(i) is not set. step 1 609 Is jumped to. while when the cumulative addi- 
45 tion request flag XCDLN(i) is set. step 6102 is proceeded to. At step 1602, the cumulative addition request flag 
XCDLN(i) is reset. Next, at step 1603, the amount of fluctuation of the torque DLN(i} is added to the cumulative value 
DLNl(i) of tiie amount of fluctuation of the torqua Next, at step 1604. the cumulative count value CDLNI(i) is incre- 
mented by exactiy 1 . That is, for example, If the cumulative addition request flag XCDLN(I) is set for the No. 1 cylinder 
at step 1601 . at step 1602, tiiis flag XCDLN(I) is reset at st^ 1603, the amount of fluctuation of the torque cumulative 
50 value DLNI(1) is calculated, and, at step 1604, the cumulative count value COLNI (1) is incremented by exactiy 1 . 

Next at step 1605, whether the cumulative count value CDLJ^JI(i) has become "8" or not is judged. When CDLNI(I) 
Is not "8'*, step 1609 is jumped to, while when CDLNI(i) becomes "8", step 1606 is proceeded to, where the torque fluc- 
tuation value DLNISM(i) of each cylinder is calculated from the following equation: 

55 DLNlSM(i) = DLNISM (i)+{DLNI (i)-DLNISM (!)}/4 

Next, at step 1607, the cumulative value DLNI(i) of tiie amount of fluctuation of the torque for each cylinder is 
deared. then, at step 1608, the cumulative count value CDLNI(i) is deared. 

That is. when there is a difference between tiie calculated amount of fluctuation of tiie torque cumulative value 
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DLNl(i) and the previously used torque fluctuation value DLNISM(0. the value of this difference {DLNI(I)-DLNISM (i)} 
multiplied by 1/4 is added to the torque fluctuation value DLNISM(i). Therefore, for example, when the cumulative count 
value CDLNI(1) becomes "8** for the No. 1 cylinder #1 , at step 1606. the torque fluctuation value DLNISM(1) is calcu- 
lated. 

5 Next, at step 1 609. whether the count value CDLNIX calculated by the routine shown in Rg. 20 has become "8" or 

not is judged. When CDLNIX is not "8", the processing cyde is ended, while when CDLNIX becomes "8". step 1610 is 
proceeded to^ where whether or not the average value of the engine speed Nave fourxi at step 1509 of the routine for 
processing CDLNIX shown In the Fig. 20 Is lower than even the set speed NL or not Is judged. When N/^ve < 
1611 is proceeded to. where the mean value of the torque fluctuation value DLNISM(i) of each torque, that is. the mean 

10 torque fluctuation value DLNISM (= {DLNISM(1)+DLNISM(2)+DLNISM(3)+DLNISM(4}}/4) is calculated. Next, step 
1615 is proceeded to. 

As opposed to this, when it is judged at step 1610 that Nave ^ NL, step 1612 is proceeded to. where whether the 
average value of the engine speed Nave '^^wer than the set speed NH (> NL) or not is judged. When Nave < NH, step 
1 61 3 is proceeded to, where tfie mean value of the torque fluctuation value DLNISM (i) of the No. 3 cylinder #3 arxi the 

IS No. 4 cylinder #4. that is. the mean torque fluctuation value DLNISM (= {DLNISM(3)+DLNISM(4)}/2) is calculated. 
Next, step 1615 is proceeded to. That is. when NL ^ Nave < NH. calculation of the amount of fluctuation of the torque 
DLN(i) for the No. 1 cylinder #1 and the No. 2 cylinder #2 is prohibited, so at this time the mean torque fluctuation value 
DLNISM is made the mean value of DLNISM(3) and DLNISM(4). 

On the other hand, when it is judged at step 1612 that Nave ^ NH. step 1614 is proceeded to. where the torque 

20 fluctuation value DLNISM (4) of the No. 4 cylinder #4 is made the mean torque fluctuation value DLNISM. That is, when 
Nave ^ NH. calculation of the amount of fluctuation of the torque DLN(i) for the No. 1 cylinder #1 . Na 2 cylinder #2. and 
No. 3 cylinder #3 is prohibited, so at this time DLNISM(4) is made the mean torque fluctuation value DLNiSM(4). Next, 
step 1615 is proceeded to. At step 1611, the count value CDLNIX is cleared. The value DUMISM representing the 
amount of fluctuation of the torque of the engine is calculated in this way. 

2S Next, the routine for calculation of the fuel injection time will be explained referring to Rg. 39. 

Referring to Rg. 39, first, at step 1801, the basic fuel injection time TP is calculated from the map shown in Rg. 2. 
Next, at step 1802, whether the operating state is one where a lean operation should be performed or not is judged. 
When tiie operating state is one where a lean operation should be performed, step 1803 is proceeded to, where the 
value of the stoichiometric air-fuel ratio feedback correction coefficient FAF is fixed at 1 .0. Next, at step 1 804, the lean 

30 con-ection coefficient FLEAN is calculated from the map shown in Rg. 4. then at step 1 805, the lean limit feedback cor- 
rection coefficient FLLFB is read from the map shown in Rg. 5. Next at step 1806, the intercylinder correction coefH- 
: cient KGTP(i) is calculated. Next, at step 1 81 1 . the fuel injection time TAU is calculated based on the following equation: 

TAU = TP - FLEAN • FLLFB • FAF • KGTP(i)+TAUV 

35 

As opposed to this, when it is judged at step 1802 that the operating state is not one in which a lean operation 
should be performed, that is. when the air-fuel ratio should be made the stoichiometric air-fuel ratio, step 1807 is pro- 
ceeded to. where tiie lean correction coefficient FLEAN Is fixed at 1 .0. Next, at step 1808, the lean limit feedback cor- 
rection coefficient FLLFB is fixed at 1.0. Next at step 1809, the intercylinder conrection coefficient KQTP(i) is fixed at 

40 1 .0. Next, at step 1810. the stoichiometric air-fuel ratio feedback correction coefficient FAF is controlled based on the 
output signal of the air-fuel ratio sensor 1 7 so that the air-fuel ratio becomes the stoichiometric air-fuel ratio. Next, step 
181 1 is proceeded to. where the fuel injection time TAU is calculated . 

Next, an explanation will be made of a third emt)odiment which enables even more accurate detection of the drive 
force generated at each cylinder or the torque generated at each cylinder. 

45 That is, the change of the speed of reciprocal motion of the piston becomes smallest near top dead center and 
becomes largest near 90"* after top dead center. In this case, if the change of the speed of reciprocal motion of the pis- 
ton becomes large, the deceleration action caused by the inertia on the piston and other reciprocating members (here- 
inafter referred to as sinply tiie inertia of tiie piston) acts on the crankshaft and therefore the change of the speed of 
reciprocal motion of the piston becomes larger near 90° after top dead center and tiie angular velocity of the crankshaft 

50 is reduced by the deceleration action due to the inertia of the piston. In this case, the higher the engine speed becomes, 
the larger the change of the speed of rec9>rocal motkm of the piston and therefore tiie deceleration action due to the 
inertia of the piston becomes more powerful the higher the engine speed becomes. 

Note, however, that when the deceleration action due to the inertia of the piston works in this way. the angular veloc- 
ity of the crankshaft near 90** after top dead center, that is the second angular velocity <db, is reduced. That is. the 

55 detected second angular velocity cob becomes the amount of Increase of the angular velocity caused by tiie combustion 
pressure plus the amount of reduction of the angular velocity due to the inertia of the piston. Therefore, to find more 
accurately the drive force generated by each cylinder or the torque generated by each cylinder, it is necessary to con- 
sider the an%>unt of change of the second angular velocity cob caused by the inertia of the pistons. This will be explained 
next refening to Rg. 40. Rg. 41 A. and Rg. 41 B. 
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As explained earlier, when the engine speed becomes higher, the deceleration action due to the Inertia of the piston 
becomes more powerful. As a result, the angular velocity of the crankshaft near 90"* after top dead center, that is, the 
second angular velocity cob, is decreased, tn other words, the elapsed time T30 near 90"* after top dead center becomes 
longer due to the deceleration action due to the inertia of the piston. Fig. 40 shows the change in the elapsed time T30 
5 at the time of high engine speed, from Rg. 40. it is learned that the elapsed times Tb(1) and Tb(3) near 90** after top 
dead center become longer. 

On the other hand, when conrtjustion is performed, the oomt)ustion pressure causes the elapsed time T30 to be 
reduced from Ta(i) to TbO). That is, when combustion Is performed, the combustion pressure causes the angular veloc- 
ity of the crankshaft to rise from the first angular velocity ooa to the second angular velocity cob. Note, however, that as 

10 explained above, the deceleration action due to the Inertia of the piston causes the second angular velocity cob to be 
decrecised, therefore the detected second angular velodty cob includes not only the amount of increase of the angular 
velocity due to the combustion pressure, but also the amount of decrease of the angular velocity due to the inertia of 
the piston. Therefore, in this case, in so far as the second angular velocity cob is not Inaeased by the amount of 
decrease of tiie angular velocity due to the inertia of the piston, the difference between the first angular velocity gol and 

IS second angular velocity cob does not express the generated drive force of the cylinder and the difference between tiie 
square of the first angular velocity coa and the squcure of the second angular velocity cob does not express the generated 
torque of the cylinder. 

Therefore, in the third embodiment, when the supply of fuel is stopped at the time of a deceleration operation, the 
difference between the elapsed times Ta(i) and TbO) >s found arKi then the difference is used to conrect the elapsed time 

20 Tb(i). That is, when the supply of fuel is stopped, no combustion pressure occurs, therefore, at this time, only the effect 
due to the inertia of the piston appears in the elapsed time T30. Rg. 41 A shows the change in the elapsed time T30 
when the supply of fuel is stopped. As shown in Fig. 41 A. when the supply of fuel is stopped, the deceleration action 
due to the inertia of the piston causes the elapsed times Tb(1) arKi Tb(3). that is. Tb(i) to become longer. At this time, 
the difference r (r= Tb(i)-Ta(i)) between the elapsed times Ta(0 and Tb(i) expresses the amount of increase of the 

25 elapsed time Tb{0 due to the inertia of the piston. Therefore, by subtracting this amount of increase r from the actually 
detected Tl)(i) shown in Rg. 40, tills result of subtraction (actually detected Tb(i)'amount of increase r) comes to 
express the elapsed time Tb(i) reduced by tiie combustion pressure. The elapsed time Tb(i) decreased by this combus- 
tion pressure is shown in Fig. 41 B. Therefore, by calculating the generated drive force arxi the generated torque of a 
cylinder based on the elapsed time Ta(i). Tb(i) shown in Rg. 41 B. the generated drive force and the generated torque 
'SO come to express the true generated drive force and generated torque from which the effect of the inertia of the pistons 
has t>een eliminated. - ' . 

In other words, at the time when the supply of fuel is stopped, the second angular velocity cob is reduced by the 
inertia of the piston t)y exactiy the difference between the first angular velocity ooa and the second angular velodty db. 
By adding the amount of decrease of the angular velocity (aa-cob) at this time to the actually detected second angular 

35 velocity cob. tiie result of the addition comes to express the second angular velocity cob Increased by the conlbustion 
pressure. 

Note, however, that the amourrt of increase r of the elapsed time due to the above inertia of the piston becomes a 
function of the engine speed. That is. if the engine speed becomes high, the deceleration action due to the inertia of the 
piston becomes powerful, so the amount of increase r of the elapsed time becomes larger, but the elapsed times Ta(i) 
40 and TbO) themselves become shorter, so the amount of increase r of the elapsed time becomes smaller as the engine 
speed N becomes higher as shown by the broken line in Rg. 42. Therefore, in this emtxKliment according to the present 
invention, the amount of increase r of the elapsed time in accordance with the engine speed N is calculated from the 
relationship shown by the broken line in Rg. 42. The amount of increase r is used to correct the elapsed time Tb(i). 

Note, however, that a value found in advance by experiments may be used as the amount of increase r of the 
45 elapsed time, but in tiiis embodiment according to the present invention, this amount of increase r is learned. That is. 
as shown in Rg. 42, a plurality of leaming regions , S3. S4. S5 ... are set for the engine speed N. When the supply 
of fuel is stopped at the time of a decela'ation operation, tiie amount of increase r^O)* ■'aCO* i^sO)* r4(i). r5(i) ... of the 
elapsed time according to the engine speed N at that time Is found for each cylinder. 

On the other hand, in the third embodiment, the rotor 13 has an outer tooth formed every 10** crank angle on its 
so outer periphery and has some of the outer teeth removed for detection of the top dead center of compression of the No. 
1 cylinder for example. The portion where the outer teetii have been removed, that is. this non-tooth portion, is shown 
by the reference numeral 18 in Rg. 43. Furtiier, in the third emtxxJiment, tiie positional relationship between the crank 
angle sensor 14 and the non-tootti portion 18 at ATDC 60"* to ATDC 90° of the No. 2 cylinder #2 and the No. 3 cylinder 
#3 is set so tiiat the crank angle sensor 1 4 faces ttie non-tootii portion 1 8. However, if the crank angle sensor 1 4 faces 
55 the non-tooth portion 18 at ATDC 60** to ATDC 90** of the is No. 2 cylinder #2 and the No. 3 cylinder #3 in this way, the 
elapsed time Tb(2) of the No. 2 cylinder #2 and tfie elapsed time Tb(3) of tiie No. 3 cylinder #3 cannot be accurately 
detected, tiierefbre tiie amount of inaease r(2) d the elapsed time for the No. 2 cylinder #2 and the elapsed time r(3) 
for the Na 3 cylinder #3 can no longer be accurately detected. 

That is, at the non-tootii portion 18, the output signal of the crank angle sensor 14 fluctuates by a large €unount 
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according to the engine speed N. As a result, as shown in Fig. 44. the amount of increase r(i) of the elapsed time for 
the Na 2 cylinder #2 and the No. 3 cylinder #3 largely deviate from the amount of increase r(i) of the Na 1 cylinder #1 
and the No. 4 cylinder #4 showing the true values. Therefore, in the third emlxxiiment. the drive force and the torque 
generated by the cylinder cannot be accurately detected for the No. 2 cylinder #2 and the No. 3 cylinder #3. therefore 

5 the drive force and the torque generated by the cylinder are made not to be calculated for the No. 2 cylinder #2 and the 
No. 3 cylinder #3. Of course in this case, as in the first embodiment and the second embodiment, if the positional rela- 
tionship between the crank angle sensor 14 and the non-tooth portion 18 is set so that the aank angle sensor 14 is 
made not to face the non-tooth portion 18 when detecting the elapsed times Ta(i) and Tb(0. it becomes possible to 
accurately calculate the drive force and the torque for all cylinders. 

10 Note that in tiie same way as in the second embodiment, in tiie third embodiment as well, the fuel injection time 
TAU is calculated based on tiie following equation: 

TAU = TP • FLEAN • FLLFB • FAF • KQTP(0+TAUV 

IS Next the routine for finding the torque generated by each cylinder will be explained refening to Rg. 45 to Fig. 53. 
Note that. Rg. 53 shows the timing of calculation of the various values p^ibrmed in each routine. 

Rg. 45 shows an interruption routine performed at every 30* aank angle. Referring to Rg. 45, first, tiie routine (step 
21 00) for calculating tiie elapsed times Ta(i) and Tb(i) is proceeded to. This routine is shown in Rg.46. Next, the routine 
(step 2200) for checking whether the calculation of the torque Is permitted or not is proceeded to. This routine is 
20 explained previously in Rg. 13 to Rg. 15. Next the routine for calculating the torque (step 2300) is proceeded ta This 
routine is shown in Rg. 47 to Rg. 49. Next, the routine for calculating the ratios KTa(i) and KTb(i) (step 2400) is pro- 
ceeded to. This routine is shown in Rg, 50 and Rg. 51 . Next, the routine for processing of tiie counter CDLNIX used for 
calculation of the torque f lLK:tuation value (step 2500) is proceeded to. This routine is shown in Fig. 52. 

Referring to Fig. 46 showing the routine for calculation of the elapsed times Ta(i) and Tb(i), first, at step 2101. the 
25 time is made the TIMEO. The electronic control unit 20 is provided witii a free run counter for showing the tima The time 
is calculated from the count value of this free run counter. Next, at step 2102. tiie current time is fetched. Therefore, the 
TIMEO at step 101 shows the time at 30" crank angle before. 

Next at step 2103. whether the No. 1 cylinder is currentiy at ATDC 30" or not is judged .When the Na 1 cylinder is 
not currentiy at ATDC 30". step 2107 is jumped to, where whether tiie Na 1 cylinder is currentiy at ATDC 90" or not is 
30 judged. When the No. 1 cylinder is not currentiy at ATDC 90", step 2109 is jumped to. 

. As opposed to this, when it is judged at step 2103 that the No. 1 cylinder is currently at ATDC 30". step 1204 is 
-• r proceeded ta where tiie tiie final elapsed time Ta(i) from TDC to ATDC 30" of the Na 1 cylinder is cateulated based on 
• the fbllowing equation: 

35 Tafi) = KTa(i)*CTIME-TIMEO) 

That is, for example, when the No. 1 cylinder #1 is currentiy at ATDC 30". the final elapsed time Ta(1) from TDC to 
ATDC 30" of tiie Na 1 cylinder #1 is calculated from KTa(1) *(T1ME-TIME0). Here. (TIME-TIMEO) expresses tiie 
elapsed time Ta(1 ) actually measured from tiie crank angle sensor 1 4 and KTa(l ) is the ratio for correcting the error due 

40 to tiie spaces between the outer teeth of the rotor 13. therefore the final elapsed time Ta(1) ot)tained by multiplying 
(TIME-TIMEO) with KTa(l) accurately expresses the elapsed time when the crankshaft rotates by 30" crank angle. 

Next at step 2105. whether the No. (i-1) cylinder where combustion had been performed one time before is the No. 
1 cylinder #1 or the Na 4 cylinder #4 or not is judged. When the No. (i-1) cylinder where combustion had been per- 
formed one time before is the No. 2 cylinder #2 or the No. 3 cylinder #3, st^ 2109 is jumped to. As opposed to this. 

45 when tiie No. (i-1) cylinder where oontestion had been performed one time before is tiie No. 1 cylinder #1 or the No. 4 
cylinder #4. step 2106 is proceeded to. 

At step 2106, the flag XCAL (i-1) of the No. (i-1) cylinder where combustion had been performed one time before 
showing that the generated torque should be calculated is set (XCAL (i-1) ^ T). In this emtxxiiment according to the 
present invention, as explained above, the ignition sequence is 1-3-4-2. so when the No. 3 cylinder #3 is currentiy at 

so ATDC 30". the flag XCAL^I) of the No. 1 cylinder #1 where combustion had been performed one time earlier showing 
that the generated torque should be calculated is set. In the same way, as shown in Rg. 53, when the final elapsed time 
Ta(2) is calculated, tiie flag XCAL(4) is set. That is, tiie flag showing that the generated torcpje shouM be calculated is 
set only for the No. 1 cylinder #1 and the No. 4 cylinder #4. 

On the other hand, when it is judged at step 21 07 that the No. 1 cylinder is currently at ATDC 90". step 21 08 is pro- 

55 ceeded to, where the final elapsed time Tb(i) from ATDC 60" to ATDC 90" of the No. 1 cylinder is calculated from the 
following equation: 

TbO) = KTb(i) • CTIME-TIMEO) 
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That Is, for example, when the No. 1 cylinder #1 Is currently at ATDC 90**, the final elapsed time Tb(1) from ATDC 
60*" to ATDC 90** of the No. 1 cylinder #1 is calculated from KTb(1) • (TIME-TIMEO). In this case as well, since the ratio 
KTb(1) for correction of the error due to the spaces between the outer teeth of the rotor 13 is multiplied with (TIME- 
TIMEO), the final elapsed time Tb(1) accurately expresses the elapsed time when the crankshaft rotates by 30*" crank 
5 angle. 

Next, at step 2109, whether the No. 1 cylinder is currently at ATDC 210** or not is judged. When the No. 1 cylinder 
is currently at ATDC 210"*. the elapsed time T18O(0 from ATDC 30** to ATDC 210*" of the No. 1 cylinder, that is. the 
elapsed time T1 80(i) of the expansion stroke of the No. 1 cylinder is calculated. While not explained In detail in the spec- 
ification of this application, the difference between the elapsed time T1 80(1) of the expansion stroke of each cylinder and 

10 the elapsed time T1 80(0 of the next expansion stroke of the cylinder is calculated for each cylind^ and tiie intercylinder 
correction coefficient KQTP(i) is calculated for each cylinder so that the difference for each cylinder becomes smaller. 
Note that this intercylinder correction coefficient KGTP(i) is disclosed in detail in for example, Japanese Unexamined 
Patent Publteation (Kokai) No. 4-370346 as mentioned above. 

Next, the torque calculation routine shown in Fig. 47 to Rg. 48 will be explained. Referring to Fig. 47 to Fig. 48. first. 

IS at step 2301 . whether the supply of the fuel has been stopped or not at the time of a deceleration operation, that is 
whether the fuel has been cut or not is judged. When the supply of fuel is not stopped, step 2313 is jumped to. while 
when the supply of fuel is stopped, step 2302 is proceeded to. At step 2302. whether the amplitude AMP calculated in 
the routine for checking permission for calculation of the torque shown in Fig. 13 to Rg. 18 is larger than the setting Bq 
or not is judged. When AMP > Bq. step 2313 is junrped to. while when AMP ^ Bq. step 2303 is proceeded to. At step 

20 2303. whether the No. 1 cylinder is cun^entiy at ATDC 90° or not is judged. When tiie No. 1 cylinder is not cunrentiy at 
ATDC 90*. step 2313 is jumped to. As opp(^ed to this, when the Na 1 cylinder is currentiy at ATDC 90"*, step 2304 is 
proceeded to. 

At step 2304, whether tiie No. (i-1 ) cylinder which was in the expansion stroke one time before (at this time, no com- 
txjstion is performed) is the No. 1 cylinder #1 or the No. 4 cylinder #4 or not is judged. When the Na (i-1 ) cylinder which 

2S was in the expansion stroke one time before is the No. 2 cylinder #2 or the No. 3 cylinder #3. step 2313 is Jumped to. 
As opposed to this, vyhen the No. (i-1) cylinder which was in the expansion stroke one time t>e1bre is the No. 1 (finder 
#1 or the No. 4 cylinder #4, step 2305 is proceeded to 

At step 2305. whether the No. (1-1) cylinder which was in the expansion stroke one time t)efore was at one of tiie 
calculation regions . S2. S3. S4. S5 ... shown in Rg. 42 for tiie engine speed N is judged. When the No. (i-1) cylinder 

30 was not at one of the calculation regions Si, S2. S3. S4. S5 .... step 2313 is jumped to. while when the No. (i-1) cylinder 
was at one of the calculation regions St, S2. S3 . S4. S5 .... step 2306 is proceeded to. At step 2306. the cfifference 
T(i-1) (= Tb(i-l)-TaO-l)) between the elapsed time Tb(i-1) of the No. (i-1) cylinder which was in the expansion stroke 
one time before and the elapsed time Ta(i-1) Is calculated. That is. when the cylinder which was in the expanston stroke 
one time before is the Na 1 cylinder #1 . the difference T(1 ) (= Tb(1 )-Ta (1 )) is calculated, while when the cylinder which 

3S was in the expansion sti'oke one time before is the No. 4 cylinder #4, the difference T(4) (= Tb(4)-Ta(4)) is calculated. 
Next, at step 2307, the difference T(i-1) is added to the cumulative value 2:T(i-l) of the cfifference 
T(i-1) (£T(i-1) = £T(i-1)+T(i-1)) . Next, at step 2308, the count value Cn(i-I) for tfie calculation region SN of the No. (i- 
1) cylinder is incremented by exactly 1. Next, at step 2309. whether the engine speed N has passed tiie calculation 
region Sn or not is judged. When the engine speed N is in the calculation region Sn. step 231 3 is junrped to, while when 

40 engine speed N has passed the calculation region Sn. step 231 0 is proceeded ta At step 231 0. the amount of increase 
rn(i-1} of tiie elapsed time is calculated based on the following equation. 

r„(i-1) - r„(i-1)+{2:T{i-1)/C„(i-1)-r„(i.1)}/16 

45 That is. for example, when the No. (i-1) cylinder which was in the expansion stroke one time before is the Na 1 cyl- 
inder #1 and the calculation region SN is Si in Rg. 42, the amount of increase ri(1) of the elapsed time is represented 
as shown in the folfowing equation: 

ri(1) = ri(1M2:T(1)/C,(1)-ri (1)}/16 

so 

That is. £T(1)/Ci(1) shows the mean value of tiie difference T(i-1) of the elapsed time when the No. 1 cylinder #1 
is in the calculation regfon S-t . When there is a difference between this mean value £T(1)/Ci(1) and the previously used 
amount of increase r<|(1). the amount of Increase r ^(1) is updated so as to approach the mean value LT(1)/Ci(1). At 
step 310, the amount of increase rnO*1) is updated, then step 231 1 is proceeded to. where the cumulative value £T(i- 
55 1) is cleared, then at step 312, the count value Cn(i-1) 's cleared. 

On the other hand, at step 2313. whether the flag XCAL(i-1) of tiie No. {\-^) cylinder where combustion had been 
performed one time kiefore showing that the generated torque should be calculated is set or not is judged. When the 
flag XCALO-1) » "0". that is. ttie flag XCAL^i-l) is not set. the processing cyde is ended. As opposed to tiiis, when ttie 
flag XCAL(i-1 ) = "1" . that is. the flag XCAL(i-1 ) is set. step 231 4 is proceeded to. where the flag XCAL(i-1} is reset, then 
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step 2315 is proceeded ta 

At step 231 5, whether the prohibition flag XNOCAL of the cylinder at which the combustion was performed one time 
before prohibiting the calculation of the torque is reset (XNOCAL = "0*0 or not Is judged. When this prohibition flag is set 
(XNOCAL = "1"), step 2324 is proceeded to. where the prohibition flag XNOCAL is reset. As opposed to this, when the 
5 prohibition flag is reset step 231 6 is proceeded to. That is, when the flag XCAL is set and the prohibition flag XNOCAL 
is reset, step 2316 is proceeded to. 

At step 231 6, the amount of increase rn(i-l) according to the current engine speed N is calculated by interpolation 
from the amounts of increase r^O). ^zQ)* CaO). UQ)* CsO) found for the calculation regions S^, S2. S3, S4. S5 ... Next, at 
step 231 7, the amount of fluctuation h of the elapsed time due to the torsional vforation of the engine drive system (Rg. 
10 8) is calculated based on the following equation: 

h = {Ta(i-1)-Ta(i)} • 60/180 

That is, as will be understood from Rg. 8, the amount of fluctuation h of the elapsed time becomes one-third of 
15 h ^ (s Ta(i-1 )-Ta (i)) . Next, at step 23 1 8. TbXi-1 ) expressing only the elapsed time reduced due to tiie combustion pres- 
sure is calculated based on the following equation: 

TbTi-1)=Tb(l-1)-r„(l-1)+h 

20 TTiat is. when finding Tb*(1) for the No. 1 cylinder #1 . h = {Ta(1)-Ta (3)} - 60/1 80 and TbXl) = Tb(1)-r „(1 )+h . Fur- 
then when finding TbX4) for the Na 4 cylinder #4. h = {Ta(4)-Ta(2)} • 60/180 and Tb'(4) = Tb(4)-rn(4)+h . 

Next, at step 2319, the generated torque DN(i-l) of the cylinder at which the combustion was performed one time 
before is calculated based on the following equation: 

25 DN(i-1) = G>b^-ffla^=»(30»/Tb*(i-1))^-(30**/Ta(i-1)^ 

This generated torque DN(i-l) expresses the torque after the elimination of the effect due to the inertia of the piston, 
the effect due to the torsional vibration ol the engine drive system, and the effect due to the 
the outer teettnof the rotor 13. > - 

30 Note that, when finding the drive force GN(i-l) generated by the cylinders, this drive force GN(i-1) niay be csUcu- 
lated t>ased on tiie following equation: 

GN(I-I) = (30'»m)(i-1))-(30»/Ta(i-1)) 

35 At step 231 9. the generated torque DN(i-1 ) is calculated, then step 2320 is proceeded to, where the amount of fluc- 
tuation of the torque DLN(i-l) in one cycle of the same cylinder is calculated based on the following equation: 

DLNfi-1) = DN(i-l)j-DNCi-l) 

40 Here, DN(i-1)j expresses the generated torque of tiie same cylinder one cycle (720** aank angle) before for DN(i-l) 
Next, at ^ep 2321, whether the amount of fluctuation of the torque DLN(i-l) is positive or not is judged. When 
DLN(i-l) ^ 0, step 2323 is jumped to, where the cumulative addition request flag XCDLN(t*1) of the cylinder at which 
the combustion was performed one time before showing that the amount of fluctuation of tiie torque DLN(i*1 } should be 
cumulatively added is set (XCDLN(t-l) ^ "1"). /\s opposed to this, when DLN(i-l) < 0. step 2322 is proceeded to. where 

45 DLN(i-l) is made zero, then step 2323 is proceeded to. Note that, tiie torque of each cylinder repeatedly rises and falls. 
In this case, to find the amount of fluctuation of tiie torque, it is sufficient to cumulatively add either of ttie amount of 
increase of tiie torque or the amount of decrease of the torque. In the routine shown in Rg. 47 to Rg. 49, only tiie 
amount of decrease of tiie torque is cumulatively added. Therefore, as explained above, when DLN(i-1) < 0. DLN(i-l) 
is made zero. 

so Next, the routine for calculating tiie ratios KTaO) and KTb(0 will be explained refenring to Rg. 50 and Rg. 51 . 

Refem'ng to Rg. 50 and Rg. 51 , first, at step 2401 , whether the supply of fuel has been stopped during the decel- 
eration operation or not that is whether tiie fuel has been cut or not is judged. When the fuel has not been cut, step 
2415 is proceeded to, where the cumulative values i:Ta(i) and £Tb(i} of the elapsed times Ta(i) and Tb(0 are cleared, 
then the processing cycle is ended. As opposed to this, when the fuel is cut, step 2402 is proceeded to. where whether 

55 tiie amplitude AMP calculated by the routine for checking permission for calculation of the torque is larger than the set- 
ting Bq or not is judged. When AMP > Bq. step 2415 is proceeded to, while when AMP ^ Bq. step 2403 is proceeded to. 

At step 2403 to step 2408. KTa(i) is calculated. That is, at step 2403, the corresponding elapsed time Ta(i) for each 
cylinder is added to tiie cumulative value £Ta(i). For example, Ta{1) Is added to LTa(1) and Ta(2) is added to LTa(2). 
Next, at step 2404. whether the Ta(i) for each (finder has been cumulatively added n numt)er of times each or not is 
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judged. When not cumulatively added n number of times each, step 2409 is jumped to, while when cumulatively added 
n number of times, step 2405 is proceeded ta At step 2405, the mean value Ma(=» [LTa ( 1 )+£Ta (2)+i:Ta (SK^'^a WV^) 
of the cumulative values i:Ta(i) of the cylinders is calcutated. Next, at step 2406. the connection value a(Q {= Ma/£Ta(i)) 
for each cylinder is calculated. Next, at step 2407. the ratio KTa(i) is updated based on the fbllowring equation: 

5 

KTa(i) KTaO)+{a{0-KTa(OJ/4 

In this way. the ratios KTa{t), KTa(2). KTa(3). and KTa(4) for the cylinders are calculated. For example, when a(1) 
has become larger than the KTa{1) used up to then, one-quarter of the difference between ot(1) and KTa(1) {a(1)- 

10 KTa(1)} is added to KTa{1), so KTa(1) gradually approaches a(1). When the KTa(i) is calculated for each cylinder at step 
2407. step 2408 is proceeded to. where the cumulative value £Ta(i) for each cylinder Is cleared. 

On the other hand, at step 2409 to step 2415, the KTbffl for the No. 1 cylinder #1 and the No. 4 cyfinder #4 are 
calculated. That is. at step 2409. whether the Tb(1) of the No. 1 cylinder #1 or the Tb(4) of the No. 4 cylinder #4 has 
been calculated or not is judged. If the Ib^) of the No. 1 cyfinder #1 or the Tb(4) of the Na 4 cylinder #4 have not been 

75 calculated, the processing cycle is ended. As opposed to this, when the Tb(1) of the No. 1 cylinder #1 or the Tb(4) of 
the No. 4 cylinder #4 has been calculated, step 2410 is proceeded to. 

At step 241 0. the corresponding elapsed time Tb(i) is added to the cumulative value i:Tb(i) for the No. 1 cylinder #1 
or the No. 4 cyfinder #4. That is. Tb(1) is added to LTb(l). while Tb(4) is added to £Tb(4). Next, at step 241 1 . whether 
Tb(l) has been cumulatively added n number of times or not is judged. When not cumulatively added n number of times 

20 each, the processing cyde is ended, while if cumulatively added n number of times, step 2412 is proceeded to. At step 
2412. the mean value Mb(= {£Tb(1)+LTb(4)}/2) of each cumulative value £Tb(i) is calculated. Next, at step 2413. the 
confection value p(i) Mb/LTb(0) for each of the No. 1 cylinder #1 and No. 4 cylinder #4 is calculated. Next, at step 
241 4. the ratio KTbQ) is updated based on the following equation: 

25 ICTb(i) KTb(i)+{pa)-KTb(i)}/4 

In this way. the ratios KTb(l) and KTb(4) are calculated for the Na 1 cylinder #1 and the No. 4 cylinder #4. For 
example. If p(1) becomes larger than the KTb(1) used up to now, the 1/4 of the difference between p(1) and KTb(1) 
{p(1)-KTb(1)l/4 is added to KTb{1) and therefore KTb(1) gr^ually approaches p(1). When the KTb(i) is calculated at 
30 . step 241 4, step 241 5 is proceeded to, where the curnulative value LTb(i) Is cle^^ 

Next, the processing of the counter COLNIX will be explained refening to Fig. 52. 

Referring to Fig. 52, first, whether the No. 3 cylirKler #3 is cunrently at ATDC 30"" or not is judged. When the No. 3 
cylinder #3 is not currently at ATDC 30"*. the processing cycle is ended, while when the No. 3 cylinder #3 is currently at 
ATDC 30''. step 2502 is proceeded to. At step 2502. whether the conditions for calculating the torque fluctuation value 

35 Stand or not Is jiKJged. For example, when the conditions for making the air-fuel ratio lean do not stand or the amount 
of change per unit time APM of the at}solute pressure in the surge tank 3 Is more than a setting or the amount of change 
per unit time AN of the engine speed is more than a setting, it Is judged that the conditions for calculating the fluctuation 
value do not stand, while at other times, it is judged that the conditions for calculating the fluctuation value stand. 
When it is judged at step 2502 that the conditions for cateulating the f luctucition value stand, step 2508 is proceeded 

40 to, where the count value CDLNIX is incremented by exactly 1 . The increment action of this count value CDLNIX is per- 
formed each time the Na 3 cylinder #3 reaches ATDC SO"", that is, with each 720"^ crank angle. Next at step 2509. the 
average value of the engine speed Nave the average value PMave of the absolute pressure in the surge tank 3 are 
calculated In tiie period from when the Increment action of tfie count value CDLNIX is started to when the count value 
CDLNIX is cleared. 

45 On the otiier hand, when it is judged at step 2502 tiiat the corKlttions for calculating tiie fluctuation value do not 
stand, step 2503 is proceeded to, where the count value CDLNIX is cleared. Next at step 2504, the cumulative value 
DLNI(i} of the torque fluctuation value DLNO*) for each of the No. 1 cyfinder #1 and tiie Na 4 cylinder #4 is cleared, tiien 
at step 2505. tiie cumulative count value CDLNl(i) for these cylinders is cleared. 

Next, at step 2506. tiie target torque fluctuation value LVLLFB is calculated. Next, at step 2507. tiie mean torque 
so fluctuation value DLNlSM(i) of each of the No. 1 cylinder #1 and No. 4 cylinder #4 is made the target torque fluctuation 
value LVLLFB calculated from the map of Fig. 228. 

Fig. 54 shows tiie repeatedly executed main routine In tiiis main routine, first, tiie routine (step 2600) for calculation 
of the torque fluctuation value is executed. This routine is shown in Rg. 55 and Rg. 56. Next, tiie routine (step 2700) for 
calculation of tiie lean limit feedback correction coefficient FLLFB Is executed. This routine is shown in tiie previously 
55 explained Fig. 27. Next when tiie predetermined crank angle Is reached, the routine (step 2800) lor calculation of tiie 
injection time is executed. This routine is shown in Fig. 57. Next, tiie other routines (step 2900) are executed. 

Next, the routine for calculation of tiie torque fluctuation value shown in Rg. 55 and Rg. 56 will be explained. 

Referring to Rg. 55 and Rg. 56. first at step 2601 . whetiier the cumulative addition request flag XCDLN(i) showing 
that the amount of fluctuation of the torque DLN(i) shouM be cumulatively added is set (XCDLN(0 » "1 ") or not is judged. 
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When the cumulative addition request flag XCDLNO) is not set, step 2609 is jumped to, while when tiie cumulation 
request flag XCDLN(i) is set, step 2602 is proceeded to. At step 2602. the cumulative addition request flag XCDLN(i) is 
reset. Next, at step 2603. tiie amount of fluctuation of the torque DLN(i) is added to tiie cumulative value DLNI(0 of the 
amount of fluctuation of the torque. Next, at step 2604. the cumulative count value CDLNI(i) is incremented by exactiy 

5 1 . That is. for example, assuming tiiat at step 2601 tiie cumulative addition request flag XCDLN(I) is set for the No. 1 
cylinder, at step 2602, this flag XCDLN(1) is reset, at step 2603. the amount of fluctuation of the torque cumulative value 
DLNI(1) is calculated, and at step 2604 the cumulative count value CDLNI (1 } is incremented by exactiy 1 . 

Next, at step 2605. whether the cumulative count value CDLNI(i) has become "8** or not is Judged. When CDLNt(i) 
is not "8". step 2609 is jumped to, while when CDLNI(i) becomes "8". step 2606 is proceeded to. where the mean toique 

10 fluctuation value DLNISM(i) of the cylinders is calculated from tiie following equation: 

DLNISM{0 = DLNISM(i)-i-{DLNI(i)-DLNISM(0}/4 

Next, at step 2607. the cumulative value DLNI(i) of tiie amount of fluctuation of the torque for Na 1 cylinder #1 or 
75 the No. 4 cylinder #4 is cleared, then, at step 2608. the cumulative count value CDLNI(i) is reset. 

That is, if there is a difference between the calculated amount of fluctuation of the torque cumulative value DLNI(i) 
and tiie previously used amount of fluctuation of the torque DLNISM(i), the value of this difference {DLNI(i)-DLNISM(i)} 
multiplied by 1/4 is added to the amount of fluctuation of the torque DLNISM(i). Therefore, for example, when the cumu- 
lative countvalueCDLNI(l) for the No. 1 cylinder #1 t>ecomes**8".atstep2606. thetorquefluctuationvalueDLNISM(1} 
20 is calculated. 

Next, at step 2609. whether the count value CDLNIX calculated in the routine shown in Rg. 52 has become "8" or 
not is judged. When CDLNIX is not "8". tiie proc^sing cyde is ended, while when CDLNIX (becomes "8**, step 2610 is 
proceeded to, where the mean value of the torque fluctuation value DLNISM(i} of the No. 1 cylinder #1 and the No. 4 
cylinder #4. tiiat is, the mean torque fluctuation value DLNISM {DLNiSM(1)-i-DLNISM(4)}/2) . is calculated. Next, at 

25 step 261 1 . the count value CDLNIX is cleared. In this way. the value DLNISM showing the anxMint of fluctuation of the 
torque of the engine is calculated. 

Note that, as explained above, if the count value CDLNIX is increased by exactiy 1 every 720** crank angle and the 
calculation of the torque for tx>th of the No. 1 cylinder #1 and Na 4 cylinder #4 is not prohibited, when the.count value 
CDLNIX has become "8". the cumulative count values CDLNI{1) and the CDLNI(4) for tiiese cylinders will already be : 

30 Therefore, in this case, the torque fluctuation value DLNISM(i) for these cylinders is calculated. On the other hand, ' 
for example, if calculation of the torque for the Na 1 cylinder #1 is prohibited, when the count value CDLNIX has 
become "8**, the cumulative count valiie CDLNI(I) for the No. 1 cylinder #1 will not become "d**. therefore the new 
annount of fluctuation of the torque cumulative value DLNI (1) for the No. 1 cylinder #1 is not calculated. Therefore, in 
this case, when finding the mean torque fluctuation value DLNISM at step 2610. the previously calculated torque fluc- 

35 tuation value DLNISM(1) is used for just the Na 1 cylirxier#1. 

Next tiie routine for calculation of the fuel injection time will be explained referring to Rg. 57. 
Referring to Rg. 57, first, at step 2801 . tiie basic fuel injection time TP is calculated from the map shown in Rg. 2. 
Next, at step 2802, whether the operating state is one where a lean operation should be performed or not is judged. 
When tiie operating state is one where a lean operation should be performed, step 2803 is proceeded to, where the 

40 value of the stoichiometric air-fuel ratio feedback correction coefficient FAF is fixed at 1 .0. Next, at step 2804, the lean 
correction coefficient FLEAN is calculated from the map shown in Fig. 4. then at step 2805. the lean limit feedback cor- 
rection coefficient FLLFB is read from the map shown in Rg. 5. Next, at step 2806. the intercylinder correction coeffi- 
cient KGTP(i) is calculated. Next at step 281 1 , tiie fuel injection time TAU is calculated based on the following equation: 

45 TAU = TP • FLEAN • FLLFB • FAF • KQTP(i)+TAUV 

As opposed to tills, when it is judged at step 2802 that the operating state is not one where a lean operation should 
be performed, that is when the air-fuel ratio should be made the stoichiometric air-fuel ratio, step 2807 is proceeded to. 
where the lean con-ection coefficient FLEAN is fixed at 1 .0, then, at step 2808, tiie lean limit feedback correction coef- 

50 f icient FLLFB is fixed at 1 .0. Next, at st^ 2809, the intercylinder correction coefficient KGTP(i) is fixed at 1 .0. Next at 
step 2810. the stoichiometric air-fuel ratio feedback correction coefficient FAF is controlled based on the output signal 
of the air-fuel ratio sensor 17 so that the air-fuel ratio becomes the stoichiometiic air-fuel raAo. Next step 281 1 is pro- 
ceeded to, where the fuel injection time TAU is calculated. 

As already explained with reference to Rg. 29, Fig. 30A, and Fig. 308. when the engine is operated, various orders 

55 of torsional vibration occur in the crankshaft Among these, particularly tiie sixth order of torsional vibration (torsional 
vibration having 60** crank angle as cyde) has a major impact on the detection of the first angular velocity oaa. Therefore, 
in the previously explained second embodiment, to prevent the drive force generated in each cylinder and the torque 
generated in each cylinder from being mistakenly detected, tiie torque was not detected for some of the cylinders. In 
the fourth embodiment to be explained below, a method separate from the second embodiment is used to accurately 



24 



EP0 799 983A2 



find the drive force generated at each cylinder or the torque generated at each cylinder. 

That is. for the generated torque to be accurately found from the first angular velocity a>a and the second angular 
velocity db, the first angular velocity coa and the second angular velocity cob must be accurately detected. However, the 
angular velocity repeatedly fluctuates both when the first angular velocity (oa should be detected and when the second 

5 angular velocity cob should be detected, therefore in the crank angle region where the generated torque can be suitably 
detected, as the value most accurately expressing each of the first angular velocity caa and the second angular velocity 
cob. a mean value of the angular velocity in a period of a certain length must be used. In this embodiment according to 
the present invention, the mean value of this angular velocity is calculated from the elapsed time required for the crank- 
shaft to rotate by a predetermined crank angle, therefore, in this embodiment according to the present invention, to fbid 

10 the accurate first angular velocity coa and second angular velocity <db, the above-mentioned predetermined crank angle 
degree must be made a crank angle degree of a certain length. 

However, when the predetermined crank angle degree tor detecting the angular velocity is made a crank angle 
degree of a certain length in this way. a separate problem sometimes occurs. 

For example, if the crank angle range for detecting an angular velocity is made long, the latter half of the angular 

IS range for detecting the first angular vefocity <aa ends up overlapping the first half of the region of reduction (Z of Fig. 
30B) of the angular velocity co due to the sixth order of torsional vibration. As a result, this sixth order of torsional vibra- 
tion causes the first angular velocity ooa to be no longer accurately detectable. To solve this problem, it is possible to 
shift the crank angle range for detection of the first angular velocity coa to the compression stroke side so that the region 
of reductton of the angular velocity o due to the sixth order of torsional vibration is not overlapped. However, if the aank 

20 angle range for detectfon of the first angular velocity coa is shifted to the compression stroke side, the first half of the 
crank angle range ends up overlapping the crank angle region unsuitable as the region for detection of the first angular 
velocity for detection of the generated torque, therefore it is impossible to detect the accurate generated torqua 

Therefore, in the fourth embodiment according to the present invention, the crank angle range for detecting the first 
angular velocity coa is set to a crank angle smaller than the crank angle range for detecting the second angular velocity 

25 cob. Explaining this with reference to the specific example used in this emtxxjiment according to the present invention, 
the crank angle range for detecting the first angular velocity coa is made the 30*" crank angle from 20** t>efore the top 
dead center of the compression stroke (hereinafter referred to as BTDC) to 10** after the top dead center of the com- 
pression stroke (hereinafter refen'ed to as ATDC). while the crank angle range for detecting the second angular velocity . 
cob is made the 50* crank angle from ATDC 50' to ATDC 100^ 

30 If crank angle range for detecting the first angular velocity coa is set to a crank angle smaller than the crank angle 
range for detecting the second angular velocity cob in this way, the crank angle range for detecting the first angular veloc- 
ity coa will not overlap the region of reduction of the angular velocity caused by the sixth order of torsional vibration and 
further will not overlap the crank angle region unsuitable as the region of detection of the first angular velocity, so the . 
second angular velocity db of course and also the first angular velocity coa can be accurately detected. 

35 Note that, even in cases other than when the sixth order of torsional vibration occurs, there are cases where it is 
necessary to set the nrtagnitude of the crank angle range for detecting the first angular velocity ooa and the magnitude 
of the aank angle range for detecting the second angular velocity cob to different magnitudes so as to accurately detect 
the first angular velocity coa and the second angular velocity cob. For example, when the timing for detection of the sec- 
ond angular velocity cob and the timing where the crank angle sensor 1 4 faces the non-tooth portion of the rotor 1 3 over- 
do lap. to accurately detect the second £ingular velocity cob, the crank angle range for detecting the second angidar velocity 
db is made a crank angle larger than the crank angle range for detecting the first angular velocity coa. 

In the fourth emtxxiiment, the rotor 13 has an outer tooth formed every 10* crank angle on its outer periphery and 
has some of the outer teeth removed for detection of the top dead center of compression of the No. 1 cylinder for exam- 
ple. The crank angle sensor 14 produces an output pulse every time the output shaft 12 rotates by 10* crank angle 

45 except at the portion where the outer teetii have been removed, that is, this non-tooth portion. The output pulse is input 
to the input port 26. 

Next, the method for calculating the drive force generated by each cylinder and the torque generated by each cyl- 
inder will be explained using a specific example. 

First, the method of calculating the drive force generated by a cylinder and the torque generated by a cylinder will 
50 be explained with reference to Fig. 58 showing the steady state operation. As explained above, the crank angle sensor 
14 generates an output pulse each time the crankshaft rotates by 10* crank angle. Further, the crank angle sensor 14 
is disposed so as to generate an output pulse at the top dead center of the compression strote (hereinafter referred to 
as TDC) of the cylinders #1 , #2, #3, and #4. Therefore, the crank angle sensor 14 generates an output pulse with each 
1 0* crank angle giving the TDC of the cylinders #1 , #2, #3, and #4. Note that the ignition sequence of the intemal com- 
55 Ixjstion engine used in the present invention is 1-3-4-2. 

The solid line of Fig. 58 shows the elapsed time when converting the time required for the crankshaft to rotate in 
each crank angle range divided by the t>roken lines to the time required for the crankshaft to rotate 30* crank angle. That 
is. in Rg. 58, Ta(0 shows the elai::^ time T30 of the 30* crank angle from BTDC 20* to ATDC 1 0* of the No. 1 cylinder. 
As opposed to this. Tb(i) shows the time for conversion of the elapsed time T50 of tiie 50* crank angle from ATDC 50* 
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to ATDC 100** of the No. 1 cylinder to the elapsed time of 30** crank angle, that Is, three-fifths of the elapsed time T50 
of 50*» crank angle. Therefore, for example, Ta(1) shows the elapsed time from BTDC 20** to ATDC 10* of the No. 1 cyl- 
inder and Tb(1) shows three-fifths of the elapsed time from ATDC SO^" to ATDC lOO** of the No. 1 cylinder. On the other 
hand, tf each crank angle degree divided by the broken tines is divided by the elapsed time, the result of the division 

5 expresses the angular velocity ox In the fourth embodiment, 30** crank angle/Ta(i) is called the first angular velocity ooa 
of the No. 1 cylinder, while SO'* crank angle/Tb(i) is called the second angular velocity cob of the No. 1 cylinder. 

Note that the change In the elapsed time shown in Rg. 58 differs somewhat by engine, therefore the crank angle 
range for detecting the first angular velocity coa and tiie crank angle range for detecting the second angular velocity db 
are determined so that in accordance with the engine (cob-coa) best expresses the drive force generated by the engine 

10 or (ffib^-coa^) best expresses the torque generated by tiie engine. Therefore, depending on the engine, the angular 
range for detecting the first angular velocity a>a nnay be from before top dead center of the compress'on stroke BTDC 
25* to ATDC 5° and the crank angle range for detecting the second angular velocity cob may be ATDC 45* to ATDC 95*. 

Rg. 59 shows in an enleirged manner the portion where tiie elapsed time Ta(i) successively calculated for each cyl- 
irKler when the engine drive system experiences torsional vik»ration decreases. As shown in Rg. 59, the elapsed time 

IS Ta(i) decreases between Ta(1) and Ta(3) by exactly tiie time ho. The decrease in tiiis time ho may be considered to be 
due to the increase in the amount of torsion due to the torsional vibration. In this case, the amount of decreeise of the 
elapsed time due to the torsional vibration between the time Ta(1) and Ta(3) may be considered to increase substan- 
tially linearly along with the elapse of time. Therefore, this amount of decrease of the elapsed time due to the torsional 
vibration is expressed by the difference between the broken line connecting Ta(1) and Ta(3) and the horizontal line 

20 passing through Ta(1). Therefore, the torsional vibration causes the elapsed time to decrease by exacUy h between 
Ta(1) and Tb(1). 

Therefore, in tiie fourth embodiment as well, h is added or subtracted to find the elapsed time Tb'(1) decreased by 
the combustion pressure. 

Rg. 60 shows the case where the space between the outer tooth of the rotor 13 showing BTDC 20* of the No. 1 
25 cylinder #1 and the outer tooth of the rotor 13 showing ATDC 10* is smaller than tiie space between other outer teeth. 
In this case, as will be understood from a comparison of Rg. 8 and Rg. 9, the elapsed time Ta(1) ends up t>ecoming 
^ smaller than the correct elapsed time for 30* crank angle. Furtiier. at this time, as will be understood from a comparison 
of Rg. 8 and Rg. 9, tiie amount of decrease of the elapsed time due to tiie torsional vibration h' becomes smaller than : 
the correct amount of deaease h, therefore the value of Tb*(1) expressing only the elapsed time decreased by the corn- 
so bustion pressure also becomes smaller than the correct value. / 

Therefore, in the fourth embodiment as well, at the time of a deceleration operation when the engine drive system : 
doesriotexperience1orsfonalvit)ration, when tiie supply of fuel has been stopped. of 
the mean value Ta(i)m of the elapsed times Ta(i) of air of the cylinders and the mean value Tb(i)ni of the elapsed time 
Ta(i) of each cylinder and the ratio KTb(i) (= Tb(i) n\/Tb(i)) of tiie elapsed times Tb(i) of all of tiie cylinders and tiie 
35 elapsed time Tb(l) of each cylinder are found, while when the fuel Is being supplied, the actoally detected elapsed time 
Ta(0 for each cylinder is multiplied by the ratio KTa(i) to find tiie final elapsed time Ta(i) for each cylinder, while the actu* 
ally detected elapsed time Tb(i) for each cylinder is multiplied by the ratio KTb(i) to find the final elapsed time Tb(i) for 
each cylinder. 

(On the other hand, as already explained with reference to Rg. 10. the Ta(i) for each cylinder fluctuates when tiie 

40 vehicle is driving over a bumpy road. In this case, as explained atx>ve. Therefore, for example, as explained above, if 
the value of ttie amount of decrease h with respect to Tb(l ) of the No. 1 cylinder #1 is found from the inclination of the 
broken line connecting Ta(1) and Ta(3) ,the value of this amount of decrease h is calculated considerably larger than 
even the actual value. As a result Tb'(1) no longer shows the correct value .therefore the drive force and the torque gen- 
erated by the cylinder can no longer be accurately detected. When the anplitude AMP becomes large, the same occurs 

45 at the cylinder giving the smallest Ta(i) . 

Further, at the cylinder where Ta(i) changes sharply from the Ta(i) of the cylinder at which the combustion was per- 
formed one time before, the value of h deviates from the actual value, tiierefore the drive force and tiie torque generated 
by the cylinder can no longer be accurately detected. Therefore, in the fourth embodiment as well, when the amplitude 
AMP is large, the drive force or the torque for the cylinder at which Ta(i) becomes maximum or minimum is not found 

50 and further the drive force or tiie torque for the cylinder where Ta(i) changed sharply witii respect to tiie Ta(i} of tiie cyl- 
IfKler where comtxistion was performed one time before is made not to be fourvl. 

Next, tiie routine for finding the torque generate at each cylinder will be explained referring to Rg. 61 to Rg. 66. 
Rg. 61 shows the interruption routine performed at BTDC 20*. ATDC 10*. ATDC 50*. and ATDC 100*. Referring to 
Rg. 61. first the routine (step 3100) for calculating tiie elapsed times Ta{i) and Tb(i) is proceeded to. This routine is 

55 shown in Rg. 62. Next, the routine (step 3200) for checking whether calculation of tiie torque is permitted or not is pro- 
ceeded to. This routine is shown in Rg. 63 to Rg. 65. Next, tiie routine for calculating tiie torque (st^ 3300) is pro- 
ceeded to. This routine is shown in Rg; 66. Next, the routine for calculating the ratios KTa(i) and KTb(i) (step 3400) is 
proceeded to. This routine is shown in the previously explained Rg. 18 and Rg. 19. Next the processing routine (step 
3500) of the counter CDLNIX used for calculation of the torque fluctuation value is proceeded ta This routine is shown 
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in the previously explained Fig. 20. 

Referring to Rg. 62 showing the routine for calculation of the elapsed times Ta(i) and Tbffl. first, at step 3101, the 
time is made the T1ME0. The electronic control unit 20 is provided with a free run counter for counting the time. The time 
is calculated from the count value of this free run counter. Next, at step 3102, the current time is fetched. Therefore, the 

5 TIMEO of step 3101 comes to express the time when the previous interruption was performed, that is. BTDC 20''. ATDC 
10". ATDC 50^ or ATDC lOO". 

Next, at step 3103. whether the No. 1 cylinder is currentiy at ATDC 10"" or not is judged .When the No. 1 cylinder is 
not cunrentiy at ATDC 10", step 3106 is junped to, where whether the No. 1 cylinder is cun-ently at ATDC lOO"" or not is 
judged .When the No. 1 cylinder is not currently at ATDC 100"*. the routine for calculation of the elapsed time sTa(i) and 

10 Tb(i) is ended. 

As opposed to this, when it is judged at step 3103 that the Na 1 cylinder is currerrtly at ATDC 10**. step 3104 is 
proceeded to. where the final elapsed time Ta(i) from BTDC 20'' to ATDC 1 0"* of the No. 1 cylinder is calculated from the 
following equation: 

IS Tafi) = KTa© • (TIME-TIMEO) 

Here, TIMEO shows tiie time at BTDC 20^. 

That is. for example, if tiie No. 1 cylinder #1 is cunrentiy at ATDC 10'. the final elapsed time Ta(1) from BTDC 20*» 
to ATDC 10" of the No. 1 cylinder #1 is calculated from KTa(1) •(TIME-TIMEO). Here, fTlME-TIMEO) expresses the 
20 elapsed time Ta(1) of 30"* crank angle actually measured by the crank angle sensor 1 4 and KTa(1) is a ratio for correc- 
tion of tfie error due to the space between the outer teeth of the rotor 1 3 therefore the final elapsed time Ta(1 ) obtained 
t>y multiplying (TIME-TIMEO) with KTa(1) accurately expresses the elapsed time when the crankshaft rotates 30" crank 
angle. 

Next, at step 3105, the flag XCAL(i-1) of the No. fi-1) cylinder where combustion had been performed one time 
25 before showing that the generated torque shouM be calculated is set PCCAL (i-1) <r- "1 "). In this embodiment according 
to the present Invention, as explained above, the ignition sequence is 1 -3-4-2. so if the No. 1 cylinder #1 currentiy is at 
ATDC 10", tfie flag XCAL(2) of the No. 2 cylinder #2 v^ere combustion was performed one time before showing that 
the generated torque should be calculated is set. In the same way, if the final elapsed time Ta(3) is calculated as shown 
in Fig. 21 . the flag XCAL(1) is set. if ttie final elapsed time Ta(4) is calculated, the flag XCAL(3) is set. and if the final 
30 elapsed time Ta{2) is calculated, the flag XCALj(4) is set. 
■ * On the otiier hand, when it is judged at step 3106 that tiie Na 1 cylinder is cunrentiy at ATDC 100"*. step 3107 is 
proceeded to. where the final ela|:»ed time Tb(0 obtained by convaiing the elapsed time of 50*"^ crank angle from the 
ATDC 50** to ATDC 1 0O** off the Na 1 cylinder to ttie elapsed time of SO*" crank angle is calculated based on the following- 
equation: 

35 

Tb(i) = 3/5 • KTb{i) • (TIME-TIMEO) 
Here, TIMEO expresses the time at ATDC SO"". 

That is. for example, when the No. 1 cylinder is currentiy at ATDC 100®. the final elapsed time Tb(l) of the Na 1 

40 cylinder #1 is calculated from 3/5 • KTb(1) • (TIME-TIMEO). Here. (TIME-TIMEO) shows the elapsed time of the 50*» 
crank angle actually measured by the crank angle sensor 14, and KTb(1) is the ratio for correction of the error due to 
the spaces between the outer teeth of the rotor 13. Therefore, the final elapsed time Tb(1) obtained by multiplying the 
(TIME-TIMEO) by 3/5 and KTb(1) accurately expresses the elapsed time while rotating by 30'' crank angta 

Next, the routine for checking permissfon for calculation of tiie torque shown in Fig. 63 to Rg. 65 will be explained. 

45 This routine is provided for prohibiting the calculation of the torque for a specific cylinder when the amplitude AMP of 
thef luctoation of Ta(i) (Fig. 10) becomes larger due to the vehicle traveling over a bumpy road. Refen^ing to Rg. 63. first, 
at step 3201 . whetiier one of the cylinders is currentiy at ATDC 1 0" or not is judged- When none of the cylinders is cur- 
rently at ATDC 10**. the processing cycle is errcled. while when one of tiie cylinders is currentiy at ATDC 10**, step 3202 
is proceeded to. The routine from step 3202 to st^ 3223 is the same as tiie routine from step 202 to step 223 shown 

so in Fig. 15. therefore the explanation of the routine from step 3202 to step 3223 will be omitted. 

Next tiie torque calculation routine shown in Rg. 66 will be explained. Referring to Rg. 66, first, at step 3301, 
whetiier the flag XCAL(i-1 ) of the No. (i-1 ) cylinder where combustion had been performed one time before showing that 
the generated torque should be calculated is set or not is judged. When ttie flag XCAL(i-1) - "O", tiiat is. when the flag 
XCAL(i-1) is not set. the processing cyde is ended. As opposed to this, when the flag XCAL(i-1) = "1". that is, the flag 

55 XCAL(i-1) is set. step 3302 is proceeded to, where tiie flag XCAIj(i-1) is reset then step 3303 is proceeded to. 

At step 3303, whetiier tiie prohibition flag XNOCAL prohibiting the calculation of tiie torque for the cylinder at which 
the combustion was performed one time before is reset (XNOCAL s "O**} or not is judged. When tills prohibition flag is 
set (XNOC/U. = "1 "). step 3310 is proceeded to, where the prohibition flag XNOCAL is reset. As opposed to this, when 
the prohibition flag is reset, step 3304 is proceeded to. That is, step 3304 is proceeded to only when tiie flag XCAL is 
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set and the prohibition flag XNCX^AL is reset. 

At step 3304. the amount of fluctuation h of the lapsed time (Rg. 59) based on the torsional vibration of the engine 
drive system is calculated based on the following equation: 

5 h = {Ta(i-1)-Ta0)l • 80/180 

That is, as will be understood from Rg. 59. the amount of fluctuation h of tiie elapsed time becomes fburniintiis of 
ho (s Ta(i-l)-TaO)) . Next, at step 305, Tb'(i-I) expressing just the elapsed time decreased due to the combustion pres- 
sure is calculated tiased on the following equation: 

10 

TbXi-l) = Tb(i-1)+h 

That is. when finding TbX1) for the No. 1 cylinder #1. h = {Ta(1)-Ta(3)} • 80/180 and Tb'(1) = Tb(1)+h . Furttier, 
when finding TbX3) for the Na 3 cylinder #3, h » {Ta(3)-Ta(4)} • 80/180 and Tb'(3) = Tb(3)-Mi . 
15 Next, at step 3306, the generated torque DN(i-1) off the q/tinder at which tiie comtMJStion was performed one time 
before is calculated based on the following equation: 

DNO-I) = ©b^-coa^ = (30*/Tb(i-1))^-(30»/Ta(i-1))^ 

20 This generated torque DN(i-l) expresses tiie torque after elimination of the effect of the torsional vit>ration of the 
engine drive system and the effect of the fluctuation in the spaces between tiie outer teeth of the rotor 1 3, therefore this 
generated torque DN(i-l) expresses the true torque generated due to the corrftHJStion pressure. 

Note that, when finding the drive force GN(i-1) generated by the cylinder, this drive force GN(i-1 ) may be calculated 
based on tiie following equation: 

25 ■ 

GNO-I) = (30*/Tb'(i-1)M30°/Ta(i-1)) 

At Step 3306, the generated torque DN(i-l) is calculated, tiien step 3307 is proceeded to, who^e the amount of fluc- 
tuation of the torque DUM(i-l) due to one cycle of the same cylinder is calculated based on the following equation: 
30 . -• • . 

DLN{i-1) = DN(i-1)j-DNCi-1) 

... Here. DN(i-1)jexpressesthegeneratedtorqueofthesamecylinderonecyde(720^crankangle)b 

Next at step 3308. whether the amount of fluctuation of tiie torque DLN(i-1) is positive or not is judged. When 
35 DLN(i-1) ^ 0. step 3310 is jumped to. where tiie cumulative addition request flag XCDLN(i-l) of the cylinder at which 
the combustion was performed one time before showing that the amount of fluctuation of the torque DLNO^I ) should be 
cumulatively added is set (XCDLN(i-l) <- "1 ^. As opposed to this, when DLN(i-l) < 0. step 3309 is proceeded to, where 
the DLNO-I) is made 0, then step 3310 is proceeded to. Note that, the torque of each cylinder repeatedly rises and falls. 
In this case, to find the amount of fluctuation of the torque, it is sufficient to cumulatively add either of the amount of 
40 increase of the torque or the amount of decrease of the torque. In the routine shown in Rg. 66. only the amount of 
decrease of the torque is cumulatively added, therefore, as explained above, when DLN(i-i) < 0, DLN(i-l) is made 0. 

Rg. 67 shows the repeatedly executed main routine. In this main routine, first, the routine (step 3600) for calculation 
of the torque fluctuation value is executed. This routine is shown In the previously explained Rg. 25 arxJ Rg. 26. Next, 
the routine (st^ 3700) for calculation of the lean limit feedback correction coeffident FLLFB is executed. This routine 
45 is shown in the previously explained Rg. 27. Next, when the predetermined crank angle is reached, the injection time 
calculation routine (step 3800) is executed. This routine is shown in the previously explained Rg. 28. Next the other 
routines (step 3900) are executed. 

While the invention has been described by reference to specific enrdxxiiments chosen for purposes of illustration, 
it should be apparent that numerous modifications could be made thereto by those skilled in the art without departing 
so from the basic concept and scope of the invention. 

Claims 

1 . A method of detection in an internal combustion engine, comprising the steps of: 

55 

setting a first crank angle revige in a crank angle region from the end of a compression stroke to the beginning 
of an expansion stroke. 

setting a secorxi ciank angle range in a crank angle region in the middle of the expansion stroke a predeter- 
mined crank angle away from the first crank angle range. 
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detecting a first angular velocity of the crankshaft in the first crank angle range. 

detecting a second angular velocity of the crankshaft in the second crank angle range, 

finding the amount of change of the angular velocity between cylinders from the difference of the first angular 

velocity of a cylinder previously performing combustion and the f rst angular velocity of a cylinder next perform* 

5 ing combustion. 

correcting tiie second angular velocity of the cylinder previoi^ly performing the comtnistion in the direction of 
decrease when the amount of change of the angular velocity between cylinders has increased, 
correcting the second angular velocity of the cylinder previously performing the combustion in the direction of 
increase when the amount of change of the angulcir velocity between cylinders has decreased, and 

10 finding the drive force generated from a corresponding cylinder based on the first angular velocity and the cor- 

rected second angular vefodty 

2. A method of detection as set forth in daim 1 , wherein a predetermined ratio of the amount of change of the angular 
velocity between cylinders is made an amount of correction of saki second angular vefodty in a lower direction arxi 
75 higher direction. 

a. A metiiod of detection as set forth in claim 2, wherein said predetermined ratio is a ratio of the crank angles of the 
first crank angle range and tiie second crank angle range with respect to the crank angle between two consecutive 
expansion strokes. 

20 

4. A method of detection as set forth in claim 1, further having the steps of 

firxiing for each cylinder a first ratio between a mean value of the first angular velocities of all of the cylinders 
and the first angular velodty of each cylinder when the supply of fuel has stopped during engine operation, 
25 finding for each cylinder a second ratio between a mean value of tiie second angular velocities of all of the cyl- 

inders and the second angular velocity of each cylinder when the supply of fuel has stopped during engine 
operation, 

correcting tiie first angular velocity of each cylirxler by the corresponding first ratio and correcting the second 
angulcir velocity of each cylinder by the corresponding second ratio when fuel is being supplied in the engine 
30 operation. 

5. A m^hod of detection as set forth in daim 1 ^ further having the step of finding a difference of squares between a 
square of said first angular velodty and a square of said corrected second angular velocity and where sakl drive 
force shows the generated torque expressed by said difforence of squares. 

35 

6. A method of detection as set forth in daim 5, further having the step of calculating the amount of torque fluctuation 
of each cylinder from the fluctuation of the torque generated at each cylinder. 

7. A metiiod of detection as set forth in daim 6, wherein the amount of torque fluctuation of each cylirtder is expressed 
40 by the difference between the torque generated when the previous combustion was performed and the torque gen- 
erated when the next comk)ustion is performed. 

8. A method of detection as set forth in daim 6. further having the steps of 

45 finding the cumulative value of the amount of torque fluctuation by cumulatively adding successively calculated 

amounts of torque f luctuertion for each cylinder exactiy a predetermined narrber of times, 
finding a representative value of the amount of torque fluctuation for each cylinder based on said cumulative 
value of the amount of torque fluctuation, and 

finding the mean amount of torque fluctuation for all cylinders, which is the mean value of the representative 

so values of the cylinders. 

9. A method of detection as set forth in daim 8, wherein when a difference occurs between a current representative 
value and the calculated cumulative value of the annount of torque fluctuation, the representative value is updated 
so as to approach tiie calculated cunmjlative value of the amount of torque fluctuation. 

ss 

10. A method of detection as set forth in daim 6. wherein the angular velocity is detected based on an output signal of 
a crank angle sensor disposed to face the teeth of a rotor connected to the crankshaft and having non-tooth por- 
tions, the generated torque is calculated for cylinders otiier than cylinders where the crank angle sensor faces the 
non-tooth portion when detecting tiie angular velodty, and tiie amount of torque fluctuation is calculated based on 
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the calculated generated torcpje. 

11 . A method of detection as set forth in claim 6, wherein an amplitude of the fluctuation of the successively found first 
angular velocities is found and, when the amplitude becomes larger than a predetermined amplitude, the f luctua- 

5 tions of the generated torques of the cylinders where the first angular velocities were largest and smallest are pro- 
hibited from being used for calculation of the amount of torque fluctuation. 

12. A method of detection as set forth in daim 6, wherein when the difference between the first angular velocity of a 
cylinder where comtxjstton had previously been performed and the first angular velocity of a cylinder where com- 

10 bustion is next performed becomes larger than a difference fcsetween the first angular velocity of the cylinder where 
conrdDusb'on is next performed and the first angular velocity of a cylinder where combustion is performed after that, 
the use of the fluctuation of the g^erated torque of the cylinder where the combustion had k>een performed second 
is prohitMted from being used for calculation of the amount of torque fluctuation. 

IS 1 3. A method of detection as set forth in claim 1 . wherein when an amplitude of a torsfonal vibration of the crankshaft 
accompanying an increase in the engine speed increases, calculation of the drive force is prohibited for the cylinder 
at the crankshaft position where the amplitude of the torsional vibratfon becomes larger. 

14. A method of detection as set forth in claim 13, wherein calculation of the drive force is prohibited for the cylinder at 
20 the crankshaft position where the amplitude of the torsional vibration becomes larger along with an increase of the 

engine speed. 

15. A method of detection as set forth in claim 1 4. where calculation of the drive force is prohibited in the order from 
the cylinder positioned at the opposite side to a flywheel attached to the crankshaft 

25 

16. A method of detection as set forth in daim 1 . further having the Steps of 

detecting a difference of angular velocity ot>tained by subtracting from the first angular velocity the second . 
angular velocity when the supply of fuel is stopped and . 
30 correcting the second angular velocity by adding said difference of angular velocity to the second angular 

velocity when the fuel Is being supplied. 

17. A method of detection as set forth in daim 16. wherein the difference of angular velocity is calculated for every pre^ ■ 
determined engine speed region when the supply of fuel is stoi^ped and the second angular velocity is corrected 

35 based on the cfifference of angular velocity detected for the engine speed region corresponding to the engine speed 
when the fuel is supplied at that time. 

18. A method of detection as set forth in daim 16, further having the step of calculating a mean value of the difference 
of angular velocity when the supply of fuel is stopped and wherein when the fuel is being supplied, the difference 

40 Of angular speed to be added to the second angular velocity is corrected so as to approach the mean value of the 
difference of angular velocity. 

19. A method of detection as set forth in claim 1 . wherein said first crank angle range is set to a crank angle range of 
a size different from the second crank angle ranga 

45 

20. A method of detectkm as set forth in daim 19, wherein said first crank angle range is narrower than said second 
crank angle range. 

21. A method of detection as set forth In daim 19, wherein said first crank angle range extends from a crank angle 
so before top dead center of the compression stroke to a crank angle after top dead center of the compression stroke. 
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Fig. 7 
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Fig. 9 
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Fig. 11 
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Fig.12 
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Fig .13 
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Fig. 14 
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